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1.  EXECUTIVE  ABSTRACT 


This  contract  period,  March  1990  to  July  1992,  was  first  concerned  with  winding  up  the 
projects  on  the  embedded  X-ray  Mask  structure  and  on  the  'quantum  lithography'  idea. 
Both  were  published  as  journal  articles  and  in  the  PhD.  Dissertation  of  N.  Maluf  (Dept,  of 
ELectrical  Engineering,  Stanford  University,  1990).  As  a  result  of  that  earlier  work  and  of 
developments  elsewhere  it  became  clear  that  the  among  the  most  critical  issues  in  achieving 
high  precision  X-ray  masks  were  those  associated  with  achieving  high  precision  in  both 
feature  size  and  feature  placement  in  electron  beam  lithography.  Most  of  the  effort  in  this 
reporting  period  was  aimed  at  achieving  precision  in  feature  size;  notably  an  attack  on  the 
problem  of  proximity  effects.  There  were  two  approaches: 

1.  A  short  term  approach  aimed  at  correcting  proximity  effects  in  existing  electron 
beam  pattern  generators  (notably  the  ETEC  MEBES  3  and  4)  for  feature  sizes  down  to  500 
nm. 


2.  A  long  term  approach  aimed  at  avoiding  proximity  effects  by  employing  low 
energy  electron  exposure  for  feature  sizes  below  500  nm. 

The  first  approach  was  successfully  executed  during  the  first  year  of  the  contract  period 
under  a  sub-contract  to  Perkin  Elmer  EBT  (now  ETEC  Systems  Inc).  Successful  correction 
for  proximity  efffects  was  demonstrated  for  features  down  to  400  nm  for  both  sparse  and 
dense  patterns  using  a  modification  of  the  GHOST  approach.  This  technique  is  now 
available  for  use  on  MEBES  equipment.  The  report  on  this  sub-contract  is  attached  as  an 
appendix. 

The  second  approach,  the  use  of  low  voltage  (<=2KV)  electron  beam  exposure,  has  also 
been  successful.  We  have  computer  modeled  and  experimentally  demonstrated  good  feature 
size  control  for  features  smaller  than  100  nm  in  both  sparse  and  dense  regions.  For  a  10% 
change  in  dose  (far  larger  than  is  normally  allowed),  the  linewidth  variation  of  was  less 
than  lOnm.  These  experiments  were  carried  out  with  a  modified  Hitachi  S-800  scanning 
electron  microscope  operating  at  2  KV.  The  electron  optics  of  this  SEM,  while  excellent  for 
an  SEM  and  just  usable  for  experiments  in  lithography,  is  far  from  optimum  as  a  reliable 
source  of  well  focused,  low  energy  electrons.  Accordingly  we  initiated  an  investigation  into 
low  voltage  electron  optics;  this  work  is  continuing  under  a  grant  from  DARPA/(5nR. 

During  the  research  period  we  also  developed  the  idea  of  using  orientation-dependent 
etched  single  crystal  silicon  as  a  useful  calibration  plate  for  electron  beam  pattern  generators 
and  so  begin  an  attack  on  the  feature  placement  problem  in  electron  beam  mask  making. 
The  crystallographic  axes  make  for  a  naturally  square  and  straight  features  that  act  to 
complement  the  laser  interferometers,  which  give  accurate  distance  calibration,  to  ensure 
improved  pattern  fidelity.  This  project  is  being  continued  under  the  current  grant  from 
DARPA/ONR. 


2.  INTRODUCTION 

Although  there  have  been  many  ingenious  ideas  for  exposing  wafers  with  sub-optical 
resolution  nearly  all  suffer  from  the  Achilles  heel  of  inadequate  mask  technology. 
Proximity  X-ray  lithography  is  no  exception.  When  X-ray  lithography  (XRL)  was  first 
invented  minimum  feature  sizes  in  integrated  circuits  were  about  7  pm  and  it  was  envisaged 
that  X-ray  lithography  would  be  introduced  for  feature  sizes  above  1pm.  For  such  features 
satisfactory  IX  electron  beam  pattern  generation  technology  existed  for  making  the  masks. 
However  by  1989  it  appeared  that  XRL  would  be  used  priniarily  for  manufacturing  circuits 
with  feature  sizes  ^low  500  nm.  As  a  result  the  there  needed  to  be  significant 
improvements  in  the  electron  beam  technology  used  for  patterning  the  masks.  The  research 
program  supponed  by  this  contract  extended  from  December  1989  to  April  1992  and  was 
concentrated  on  the  issue  of  minimizing  feature-size  errors  introduced  by  electron  scattering 
in  the  workpiece  ('proximity  effects').  There  were  two  main  aspects: 

1.  Developing  a  short  term  correction  procedure  that  could  be  implemented  on  existing 
pattern  generators  and  was  acceptable  for  feature  sizes  down  to  500nm  and  perhaps  below. 

2.  Investigating  longer  term  strategies  that  would  suffice  for  feature  sizes  down  to  l(X)nm. 

A  third  aspect  was  related  to  the  second  and  was  an  investigation  of  advanced  electron 
optical  concepts  for  new-generation  electron  beam  mask  makers.  This  third  aspect  is  the 
main  thmst  of  the  new  contract  initiated  September  1992. 


3.  SHORT  TERM  CORRECTION  FOR  ELECTRON  BEAM  PROXIMITY 
EFFECTS 

Although  many  schemes  for  correcting  proximity  effects  had  been  described  since  the 
original  work  of  Parikh  in  1977  (Int.  Symp.  on  Electron,  Ion  and  Photon  Beam 
Tehcnology,  Proceedings  published  in  J.  Vac.  Sci.  Tech.).  None  had  been  implemented  on 
the  most  widespread  mask  making  tools,  the  MEBES  series  of  electron  beam  pattern 
generators  from  ETEC.  Accordingly  ETE<2  undertook,  as  a  1-year  sub-contract,  to  review 
die  options  and  choose  and  implement  the  most  suitable.  After  the  initial  review  a  modified 
version  of  the  GHOST  (G.  Owen  and  P,  A.  Rissmann,  J.  Appl.  Phys.  June  1983)  strategy 
was  chosen.  Computer  simulation  indicated  that  by  using  GHOST  for  correcting  for  the 
effects  of  electrons  backscattered  from  the  substrate  and  feature-size  biassing  to  correct  for 
forward  scattering  in  the  resist  it  should  be  possible  to  achieve  adequate  correction  for 
feature  sizes  down  to  below  400  nm  without  changing  the  basic  mask  making  materials  and 
processes .  This  prediction  was  borne  out  by  experiment;  as  more  experience  was  gained,  it 
was  even  found  unnecessary  to  bias  the  feature  edges  (fig.  3.1).  This  procedure  has  now 
been  implemented  as  an  option  for  commercial  mask  vendors.  The  ETEC  report  describing 
this  work  is  attached  as  an  appendix. 
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4.  ELIMINATION  OF  PROXIMITY  EFFECTS-SIMULATION 


The  major  goal  for  this  work  was  to  establish  the  ability  to  predict  the  proximity 
effects  for  electron  beam  lithography  over  a  wide  range  of  energies  and  materials.  The 
processes  in  electron  scattering  are  non-linear  and  at  medium  and  high  energies  (S- 
lOOkeV)  are  dependent  on  the  three  dimensional  structure  of  the  target.  For  this  reason 
the  simulation  of  the  deposited  energy  distribution  by  solution  of  a  difiusion  equation  is 
only  possible  with  considerable  simplifying  assumptions  and  Monte-Carlo  simulations  of 
electron  trajectories  is  the  major  technique  for  determining  the  properties  of  the  electron 
beam  target  interaction. 

For  simulation  of  electron  trajectories  we  initially  used  software  packages  obtained 
from  U.C.Berkley  (SAMPLE)  and  the  University  of  Tennesse.  We  found  that  although 
these  packages  were  excellent  for  the  conventional  range  of  electron  energies  expected  in 
electron  beam  lithography  (10-30kkeV)  they  did  not  give  us  the  energy  range  or  the 
flexibilty  in  materials  we  required.  In  the  absence  of  software  packages  available  to  fit  our 
needs  we  have  developed  a  complete  software  package  to  simulate  exposure  and 
development. 

The  simulator  software  we  have  developed  is  flexible  and  accurately  covers  a  very 
wide  range  of  target  structures  and  energies.  Multiple  target  layers  are  possible  and  thus 
electron  scattering  from  structures  such  as  multi-layer  X-ray  masks  can  be  simulated.  The 
radial  distribution  of  deposited  electron  energy  is  treated  in  a  novel  way  to  give  high 
accuracy  but  minimum  computation  time  over  a  wide  range  of  energies.  In  a  conventional 
Monte-Carlo  program  the  distribution  of  electrons  scattered  into  a  resist  is  treated  by 
considering  the  energy  deposited  into  a  set  of  parralell  stripes  (bins)  linearly  increasing 
away  from  the  center  of  the  initial  irradiation.  Ten  thousand  stripes  are  needed  to  cover 
the  range  of  the  electrons  (10  microns)  with  Inm  acuracy.  This  large  number  of  stripes 
increases  the  computation  time  both  in  the  Monte-Carlo  simulation,  because  relatively  few 
electrons  reach  the  outer  stripes,  and  also  during  the  fitting  of  Gaussians  to  the  final 
energy  distributions.  As  an  alternative  to  linear  stripes  we  have  used  a  radial  bin  with  an 
exponentially  increasing  bin  size.  An  initial  circular  bin  of  Inm  diameter  is  surrounded  by 
annular  bins  increasing  by  1 0%  per  bin.  Thus  the  first  annular  bin  is  between  1  nm  and 
l  .lnm  diameter,  and  the  bins  increase  to  10.4  microns  in  only  97  steps.  This  scheme  thus 
gives  high  spatial  accuracy  where  it  is  neeeded  at  ihe  impact  point  and  to  simulate  low 
energy  scattering,  yet  gives  a  wide  range  for  high  energy  electrons  with  large  bins  where 
the  electron  energy  deposition  per  unit  area  is  low.  We  have  found  this  scheme  to  be  very 
successful!. 

The  resist  layer  which  is  normally  chosen  as  the  top  layer  is  divided  into  five  equal 
slices  and  up  to  six  Gaussians  can  be  fitted  to  the  energy  distributions  in  each  slice.  We 
normally  find  that  at  least  three  Gausssians  are  needed  to  fit  the  resist  energy  distributions 
to  better  than  1 0%  accuracy  and  that  the  forward  scattering  is  not  well  represented  by  a 
single  Gaussian.  This  latter  fact  is  obvious  from  the  form  of  the  Mott  and  Rutherford 
scattering  crossections. 
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crossections.  We  have  for  low  energy  elenrons 

"“rtaiS"  *  atomic  number  mrgets  both  for  the  to.^  and  the  differential  elastic 
cross-section.  This  work  is  fully  described  in  off 


Low  Energy  Electron  Ineleastic  Mean  Free  Paths  In  Resists. 

The  inelastic  mean  free  path  (IMFP)  of  low  energy  (1.4keV)  electrons  traveling 
in  the  Langmuir-Blodgett  (L-B)  PMMA  films  was  measured  using  x-ray  photoelectron 
spectroscopy  (XPS).  In  an  XPS  experiment,  the  probability  of  a  primary  photoelectron 
escaping  from  a  depth  z  below  the  surface  of  a  clean  solid  A  is  proportional  to 
exp(-z  /  )  where  is  the  inelastic  mean  free  path  for  the  electrons  in  the  solid[l,2]. 

Therefore,  the  intensity  ratio  of  the  substrate  XPS  signal  ( Ifj)  with  an  overlayer  A  of 
thickness  d  to  that  (  Iq)  in  the  absence  of  the  overlayer  is  . 

Ij  =  Iq  exp(-d  /  A^  sin  0) 

With  the  knowledge  of  the  precise  thickness  of  overlayer  and  the  take-off-angle 
the  IMFP  can  be  obtained  using  the  relationship  between  measured  I^j  and  Iq.  An 
empirical  equation  based  on  a  compilation  of  results  is  often  used  for  the  IMFP  of 
electron  energies  in  the  range  between  10  to  2000eV.  For  organic  overlayers  the  IMFP 


A  = -^Y-  +  0.4l(aE)^  monolayer 

where  a  is  the  monolayer  thickness  in  nm  and  E  is  the  electron  energy  in  eV[3]. 

From  this  equation  it  can  be  seen  that  XPS  with  electron  energies  in  the  range  of 
0. 1  to  2keV  is  a  surface  sensitive  technique  with  the  majority  of  the  XPS  signal  coming 
from  the  first  few  surface  monolayers.  For  IMFP  measurements  using  XPS,  it  is  clearly 
important  to  know  the  precise  overlayer  thickness.  For  organic  molecules,  spin  casting  is 
the  common  technique  to  have  uniform  films  on  substrates.  However,  as  the  thickness 
goes  thinner  than  lOnm,  the  pinhole  density  gets  much  higher.  The  Langmuir-Blodgett 
technique  is  a  method  of  preparing  ultra-thin  L-B  PMMA  films  of  a  precise  thickness  with 
very  few  pinholes[4]. 


Atactic-PMMA  obtained  from  Pressure  Chemical,  with  weight  average  molecular 
weight  (Mw )  of  188,100  amu  and  Mw/Mn  <  1.08,  where  Mn  is  the  average  molecular 
weight,  was  used  to  prepare  the  monolayers  on  the  thermally  evaporated  Au/Si  substrate 
by  the  standard  Langmuir-Blodgett  technique.  A  Joyce-Loebl  Langmuir  Trough  IV 
equipped  with  a  microbalance  for  measurement  of  the  surface  pressure  by  the  Wilhelmy 
plate  method  was  used  to  prepare  the  Langmuir-Blodgett  films.  Filtered  deionized  water 
with  a  pH  of  7  was  used  for  the  subphase.  PMMA  was  spread  on  the  water  surface  from  a 
dilute  chlorobenzene  solution  (10  mg  PMMA  in  20  ml  benzene).  Transfer  pressure  of  15 
dyn/cm  was  used  to  deposit  L-B  films  during  downstroke  and  upstroke  at  the  speed  of  2 
mm/min.  The  first  monolayer  is  transferred  to  Au/Si  substrate  by  dipping  the  substrate  to 
the  subphase  due  to  the  hydrophobic  nature  of  the  Au  surface.  The  subsequent  monolayers 
were  deposited  during  upstrokes  and  downstrokes.  The  thicknesses  of  the  L-B  films  can 
be  accurately  measured  by  ellipsometry  technique.  For  each  L-B  PMMA  monolayer,  the 
thickness  is  0.85nirr4]  so  that  in  using  overlayer  structures  from  two  to  eight  monolayers 
overlayer  thicknesses  are  from  1.7  to  6.8  nm. 

A  Surface  Science  ESCA  system  was  used  for  the  XPS  measurements. 
Monochromatic  X-rays  from  A1  target  was  used  to  minimise  damage  to  the  polymer 
structures.  C Is,  01s,  and  Au4f  were  collected  from  the  L-B  PMMA  films  of  2,  4,  6,  and  8 
monolayers  at  two  take-off-angles;  38  and  60  degrees  from  the  surface  plane.  For  these 
experiments,  the  PMMA  L-B  films  were  prepared  separately. 

The  photoelectron  (1.4keV)  intensities  of  Au4f7/2  peaks  at  different  take-off- 
angle?  are  plotted  as  functions  of  number  of  the  L-B  PMMA  monolayer,  as  shown  in 
figure  The  IMFP  can  be  obtained  by  using  the  relationship  between  and  Iq,  as 
described  in  equation  (2).  The  IMFP's  derived  from  figure  1  are  listed  in  table  1,  and  the 
average  IMFP  for  1.4keV  electrons  traveling  in  the  L-B  PMMA  films  is  found  to  be 
8.3nm. 

1.  W.  J.  Carter,et  al,  J.  Electr.  Spectr.  5,  827(1974). 

2.  T.  A.  Carlson  and  G.  E.  McDuire,  ibid,  1,  161(1972). 

3.  M.  P,  Seah  and  W.  A.  D.  Dench,  Surf  Interface  Anal.  1,  2(1979) 

4.  S.  W.  J.  Kuan,  C.  W.  Frank,  C.  C.  Fu,  D.  R.  Allee,  P.  Maccagno,  and  R.  F.  W. 

Pease,  J.  Vac.  Sci.  Technol.  B6,  2274(1988). 


Table  1  Measured  IMFPs  at  different  take-off-angles. 


4.  I 

The  photoelectron  intensity  of  Au4f7/2  peaks  at  different  take  off  angles  as  a  function 
of  number  PMMA  L-B  monolayers,  (a)  at  38^  (b)  at  60°  take  off  angle. 
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(a)  IMFP  measurniKnt  at  take-off-angle  38*. 
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(b)  IMFP  measuTement  at  take-off-angle  60*. 


5.  ELIMINATION  OF  PROXIMITY  EFFECTS-EXPERIMENT 


Initial  experiments  to  verify  the  predictions  of  the  previous  section  were  carried  out  on  a 
Nanometrics  Cwikscan  HIE  scanning  electron  microscope  (SEM).  However  the  availability 
of  a  Hitachi  S800  SEM  offered  the  possibility  of  superior  gun  performance  at  low  voltages 
and  so,  despite  the  lack  of  a  blanking  capability,  was  used  for  the  experiments  descried 
here.  Most  experiments  were  carried  out  at  2  KV  because  the  simulations  predicted  that  an 
ener^  of  2KeV  was  a  suitable  choice  as  the  maximum  energy  to  keep  the  proximity  effects 
negligible  and  the  minimum  value  to  penetrate  a  thick  enough  resist  layer  to  allow  useful 
pattern-transfer.  The  resist  thickness  was  nominally  100  nm  but  a  more  accurate 
measurement  indicated  that  the  corect  value  was  66nm.  This  is  quite  thick  enough  to  allow 
patterning  of  50nm  chromium  which  could  be  used  directly  as  an  absorber  for  deep 
ultraviolet  light  or  as  an  etch  mask  '^or  silicon  or  silicon  dioxide  (S.  W.  J.  Kuan,  PhD. 
dissertation,  Stanford  University,  1988).  The  resist  material  was  Poly(methyl  - 
methacrylate)  (PMMA)  of  Mw=188,100=1.08Mn.  The  substrates  used  were  oxidized 
silicon  (to  simulate  an  X-ray  mask  when  using  an  additive  pattern  transfer  process)  and 
oxididized  silicon  with  300nm  gold  (to  simulate  an  X-ray  mask  proces  using  subtractive 
pattern  transfer). 

The  test  pattern  used  contained  an  isolated  exposed  line,  equal  lines  and  spaces,  and  an 
isolated  unexposed  line  (fig.  5.1).  The  exposure  levels  were  varied  from  ISiiC/cm^  to 
35pC/cm2-  As  predicted  and  previously  reported  the  required  dose  is  inversely  proportional 
to  the  primary  electron  energy  thus  compensating  exactly  to  the  loss  of  beam  brightness  at 
lower  energies  for  a  given  convergence  angle. 

The  results  showed  that  the  proximity  effect  is  indeed  virtually  eliminated  under  the  above 
conditions  for  linewidths  down  to  less  than  100  nm  (fig.  5.2).  By  measuring  the  linewidths 
as  the  fullwidth  at  half  maximum  of  the  SEM  video  signal  we  also  plotted  widths  of 
exposed  isolated  lines  and  of  unexposed  isolated  lines  as  a  function  of  exposure  for  two 
different  substrates  (figs.  5.3).  From  these  results  we  determined  that  the  linewidth 
variation  as  a  result  of  a  10%  exposure  variation  (more  than  twice  that  normally  allowed)  is 
only  6nm  and  9nm  for  the  two  representative  substrates  (table  5.1).  The  higher  figure  is  for 
the  gold  substrate  as  expected  due  to  the  increased  backscattering  from  the  higher  atomic- 
number  substrate.  These  figures  are  within  the  limit  of  any  rational  specifications  for 
feature  sizes  down  to  lOOnm. 

Since  the  use  of  the  SEM  video  signal  as  a  measure  of  linewidth  is  open  to  some  ambiguity 
we  also  investigated  the  use  of  th'*  atomic  force  microscope  (AFM)  as  a  metrology  tool. 
This  is  describe  in  the  following  section. 


6  METROLOGY  WITH  THE  ATOMIC  FORCE  MICROSCOPE. 

An  atomic  force  microscope  (AFM)  has  been  used  to  study  the  latent  image  in  an 
electron  beam  resist,  polymethylmethacrylate  (PMMA).  PMMA  shrinks  under  electron 
bombardment  and  the  dose  dependence  of  this  process  can  be  characterized  using  optical 
ellipsometry.  Knowledge  of  the  dose  dependence  of  the  shrinkage  gives  the  electron 
energy  deposition  distribution  from  the  AFM  profile.  The  half  width  of  the  "writing 
electron  beam  can  also  be  estimated  from  the  AFM  profile  giving  an  unique  measurement 
method  for  this  difficult-to-measure  quantity. 


lOO^m 


two-address  (nominal  S6nm)  lines  and  spaces 


two-address  line  in  a  dense  area 
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Fig.  JS  Proximity  effect  pattern  consists  of 
two-address  lines  in  dense  and  sparse  areas. 

The  dark  areas  are  the  exposed  areas.  The  white 
areas  are  the  unexposed  areas. 
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Tabled]  Comparison  of  meas'ired  and  modeled  linewidth  difference  at 
10%  dose  variation  for  two-address  lines  in  dense  and  sparse  area  on 
both  Si  and  3(X)nm  thick  gold  coated  Si  substrates. 


lOOnm  features  were  resolved  in  60nin  PMMA  resist  both  Si  (left)  and  300nm  Au/Si  (right)  substrates. 
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Fig.  Lincwidths  patterned  on  bare  silicon  as  a  function  of  exposure  dose  for  86nni 
features.  Note  (a)  the  absence  of  any  effect  of  the  density  of  the  pattern  and  (b)  the 
low  sensitivity  of  linewidth  to  dose  (oL=40nm  for  a  2-fold  increase  in  dose),  (c)  the 
difference  of  the  linewidths  in  dense  and  sparse  areas  is  less  than  8nm. 


Fig.  Linewidths  patterned  on  gold  coated  silicon  substrates  as  a  function  of 
exposure  dose  for  nominal  86nm  features.  Note  (a)  negligible  proximity  effect  on  these 
patterns,  (b)  the  difference  of  the  linewidths  in  dense  and  sparse  areas  is  less  than 
8nm,  and  (c)  these  patterns  can  be  resolved  over  a  two-fold  exposure  range. 


Several  recent  studies  have  suggested  the  possibility  of  using  the  atomic  force 
microscope  (AFM)  as  a  metrology  tool  for  electron  beam  resists.  In  the  study  reported 
here  we  have  used  the  AFM  to  study  the  latent  image  in  a  resist,  that  is,  after  electron 
beam  exposure  but  before  the  normal  development  of  the  resist.  The  aim  of  this 
experiment  being  to  measure  the  exposed  dose  profile  in  the  resist. 

The  resist  used  was  82nm  of  polymethylmethacrylate  (PMMA)  spun  onto  a  Si 
wafer  with  a  pre-exposure  bake  of  1  hr  at  120*^C.  All  exposures  used  a  MEBES-II  electron 
beam  lithography  tool  which  has  been  modified  to  give  both  SkeV  and  the  standard  lOkeV 
operation.  The  SkeV  exposures  are  limited  to  0.25  micron  minimum  beam  size,  while  the 
minimum  lOkeV  beam  size  is  0. 12S  micron. 

A  matrix  of  line  exposures  was  run  using  three  different  doses  60,  180  and 
600uC/cm2.  The  nominal  PMMA  exposure  dose  is  bOuC/cm^.  Both  0.125  and  0.25 
micron  single  address  width  lines  were  exposed  at  lOkeV,  and  0.25  micron  lines  at  SkeV. 
An  AFM  micrograph  over  an  8x10  micron  area  of  the  0.125  micron  lines  exposed  at 
lOkeV  and  bOOuC/cm^  is  shown  in  figure  1.  The  scale  is  exagerated  in  the  vertical 
direction.  The  exposed  lines  are  approximately  lOnm  deep  which  represents  a  12% 
shrinkage  of  the  film.  As  can  been  seen  there  is  a  great  deal  of  detail  in  this  image,  both  in 
cross-section  and  in-plane.  Most  of  the  features  seen  are  'real',  not  noise  in  the  recording 
system,  and  give  a  wealth  of  information  about  the  exposure  process  and  the  film 
topography.  Some  small  dust  particles  can  be  seen,  and  some  damage  due  to  bounce  of 
the  Si02  cantilever  over  one  of  these  is  apparent  as  a  line  on  the  trailing  edge  of  the 
particle.  In  other  AFM  scans  the  recording  force  of  approximately  lOOnN  has  been 
reduced  to  lOnN  to  avoid  such  effects. 

Profiles  across  the  0.25  micron  lines  exposed  at  5  and  lOkeV  and  bOOuC/cm^  are 
shown  in  figure  2.  The  profiles  have  been  plotted  inversely,  as  a  positive  shrinkage,  to 
emphasize  that  we  intend  to  derive  exposure  dose  distributions  from  these  profiles.  From 
both  these  figures  the  effects  of  scattering  in  the  resist  and  backscattering  from  the  Si 
substrate  can  be  clearly  seen.  For  the  5keV  exposure  the  greater  width  of  the  forward 
scattering  and  the  confined  nature  of  the  backscattering  has  produced  an  overlap  of  the 
dose  distributions  that  has  effectively  broadened  the  linewidth.  For  the  lOkeV  exposure 
the  forward  and  backscattered  distributions  are  well  separated  and  the  classic  'double 
Gaussian'  distribution  can  be  observed. 

The  results  can  be  quantified  within  limits  largely  set  by  the  AFM  accuracy  and  the 
resist  surface  roughness.  The  AFM  (Park  Scientific  Instruments)  used  a  piezoelectric 
scanner  with  a  10  micron  maximum  scan.  Thus  these  results  are  at  the  limit  of  the  scan 
size  and  do  show  some  non-linear  behavior.  The  resist  shrinkage  shown  in  the  figures  here 
are  nominal  nm  units,  but  were  not  calibrated  for  either  magnitude  or  linearity  for  the 
distance  range  used  in  these  experiments.  A  further  limit  on  the  accuracy  of  all  the 
measurements  using  the  line  profiles  is  the  roughness  of  the  PMMA  surface  (2nm  RMS) 
and  the  number  of  AFM  image  line  scans  that  can  summed  without  either  introducing 


broadening  of  the  features  or  adding  noise  from  defects.  The  profiles  shown  in  figure  2  are 
summed  over  10  linescans  which  is  equivalent  to  approximately  0.4  micron  in  the  direction 
of  the  line  exposures. 

A  further  complication  in  quantification  is  the  non-linearity  and  scan  distortion  in 
the  vertical  direction  due  to  the  limited  scan  of  the  particular  AFM  used  here.  A  linear  and 
quadratic  term  have  been  removed  fi-om  the  linescans  to  give  an  approximately  flat  base 
line,  but  there  is  necessarily  some  informality  about  this  procedure,  and  for  the  wider 
lOkeV  backscattered  distributions  this  may  introduce  a  significant  error. 

To  recover  the  dose  distribution  from  the  AFM  profiles  we  have  used  a  model  of 
resist  shrinkage  based  on  ellipsometry  measurements.  For  these  measurements  2x2mm 
areas  were  exposed  with  doses  up  to  1600uC/cm^  at  lOkeV  (2S  times  the  nominal 
exposure  dose  for  PMMA)  and  up  to  SOOuC/cm^  for  5keV  exposures.  The  dependence  of 
the  resist  thickness  with  dose  and  the  change  in  the  refractive  index  were  both  measured. 
The  results  of  the  resist  thickness  measurements  are  presented  in  figure  3  as  the 
normalized  resist  shrinkage,  (to-t)/to,  where  to  is  the  initial  film  thicknes  and  t  is  the 
exposed  thickness.  The  maximum  shinkage  was  24%  for  a  dose  of  1600uC/cm^  at 
lOkeV,  while  the  refractive  index  changed  from  1.486  to  1.509  over  this  range.  The 
change  in  the  refractive  index  can  only  account  for  20-25%  of  the  resist  shrinkage  and  the 
remainder  is  probably  volatized.  This  is  also  evidenced  by  the  pressure  rise  observed  in  the 
MEBES  lithography  tool.  The  combined  effect  of  the  removal  of  material  and  the 
densification  of  the  resist  is  to  produce  a  complex  dependence  of  the  shrinkage  on  dose. 
There  are  several  distinct  processes  that  may  contribute  to  the  shrinkage^,  scissioning  of 
the  polymer  chains  followed  by  relaxation,  liberation  of  the  various  volatile  components, 
each  with  it's  own  rate,  and  possibly  recombination  at  broken  bonds.  We  might  therefore 
expect  that  the  shinkage  can  be  characterised  by  a  summation  of  the  different  reaction 
rates  removing  that  part  of  the  resist  that  is  available  for  the  interaction  with  the  electron 
beam.  The  ellipsometer  data  of  the  present  study  shown  in  figure  3  can  be  fitted  to  within 
+0.2nm  using  two  linear  rates  and  limits  for  the  5  and  lOkeV  exposures. 


1-t/to  =  0.021(1  -  e’O  012d  )  +  o.42(l  -  e-0  00065d  )  (1) 

where  t  is  the  exposed  resist  thickness,  to  the  original  resist  thickness,  and  d  is  the  dose  in 
uC/cm2.  The  pre-exponential  factors  represent  the  proportion  of  the  resist  shrinkable  by 
one  mechanism,  and  the  exponent  the  rate  of  the  process. 

The  profilometer  results  of  Erasmus^  using  0.5  micron  PMMA  show  that  the 
shrinkage  tends  to  limit  at  0.48  the  original  thickness  with  doses  above  lOOOC/m^  and  the 
results  for  lOkeV  from  Erasmus^,  which  are  also  plotted  on  figure  3,  require  additional 
slower  rate  processes  for  agreement  out  to  these  large  doses  but  are  otherwise  consistent 
with  this  work. 

As  the  ellipsometer  shrinkage  measurements  from  82nm  thick  resists  presented 
here  are  consistent  up  to  600uC/cm2  with  the  profilometer  results  of  Erasmus^  using  a  0.5 


micron  resist,  it  would  be  expected  that  there  would  be  no  non-linearities  associated  with 
diffusion  of  volatiles  through  the  film.  Thus  given  the  apparently  high  diffiisivity  of  these 
volatiles  there  seems  less  likelihood  of  a  build-up  of  material  just  out  of  the  electron  beam. 

In  the  range  of  doses  reported  here  the  shrinkage  given  by  equation  1  is  also  seen 
at  other  energies  and  for  microscopic  structures.  Figure  4  shows  linescans  for  SkeV 
exposures  at  three  different  doses.  The  peak  heights,  which  are  plotted  in  figure  3  follow 
the  same  dose  dependence  as  the  results  for  the  lOkeV  area  exposures.  This  can  be 
understood  if  the  backscattering  factor  is  near  unity  as  then  the  peak  height  of  the  SkeV 
line  exposure  will  be  at  approximately  the  same  deposited  energy  as  the  lOkeV  area 
exposure.  The  results  of  Erasmus^  show  that  the  shrinkage  scales  approximately  linearly 
with  deposited  energy  and  are  consistent  with  our  Monte  Carlo  simulation  which  gives  a 
ratio  of  1:1.85:2  for  5,  10,  and  20keV  exposures  in  0.5  micron  resists. 

The  half  width  of  the  forward  scattering  distribution  for  the  0.25  micron  lines  at 
lOkeV  from  figure  2b  is  0.26  micron,  which  equates  to  0.23  micron  using  the  dose 
dependence  of  equation  1.  The  backscattering  distribution  is  found  to  be  0.7  of  the  dose  of 
the  forward  distribution  (assuming  no  overlap  between  lines).  The  narrower  than  expected 
forward  distribution  is  in  fact  consistent  with  the  MEBES  method  of  estimating  the  beam 
profile  width  by  measuring  across  an  edge.  Apart  from  the  error  associated  with  the  non¬ 
linear  scan  of  this  particular  AFM,  this  method  represents  a  more  realistic  estimate  for  this 
notoriously  difficult-to-measure  quantity. 

The  half  width  of  the  5keV  0.25  micron  profile  is  0.49  micron  on  the  linear  plot  of 
figure  2a  which  equates  to  a  0.42  micron  dose  distribution.  The  line  broadening  observed 
experimentally  is  consistent  with  the  overlap  of  the  forward  and  backscattered 
distributions  from  Monte-Carlo  calculations. 

The  results  reported  here  demonstrate  that  direct  measurements  of  the  electron 
beam  dose  profiles  in  a  resist  can  be  made  using  an  AFM.  Besides  giving  a  measure  of  the 
electron  beam  half  width,  the  technique  can  be  used  to  quantify  proximity  effects,  and  any 
dose  distribution  on  a  free  surface  that  can  be  coated  with  PMMA. 
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AFM  micrograph  of  the  latent  image  of  e-beam  exposed  lines. 
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7.  BRIEF  REVIEW  OF  PERFORMANCE  OF  LOW  VOLTAGE  ELECTRON 
OPTICS 

There  is  now  a  substantial  body  of  work  pointing  out  the  that  the  penalty  of  moving  to  low 
energy  electrons  (increased  minimum  beam  diameter  and  reduced  available  current  density) 
is  much  less  than  indicated  in  traditional  electron  c^tical  analysis. 

There  have  been  several  developments  leading  to  this  change  of  mind: 

1.  The  possibility  of  using  a  projection  lens  in  which  the  woikpiece  is  immersed  in  the 
focusing  field. 

2.  In  particular  the  use  of  a  retarding  field  as  the  focusing  field. 

3.  The  use  of  a  retarding  field  followed  by  a  field-free  region  for  hte  workpiece 
(Meisburger  et  al.  EIPBT  symposium,  Orlando,  Fla.,  June  1992) 

4.  The  introduction  of  low-energy  -spread  electrtMi  emitters  such  as  thermal-field  emitters. 

Although  we  have  not  yet  determined  the  optimum  configuration  for  a  low  voltage  electron 
beam  column  for  lithography  it  is  clear  that  we  can,  using  some  of  the  above  developments 
obtain  some  very  attractive  results. 

For  example  if  we  assume  that  we  employ  a  thermal  field  emitter  with  an  energy  spread  dV 
of  1  eV  and  virtually  unlimited  brightness  then,  neglecting  space  charge  we  can  estimste  the 
minimum  beam  diameter  obtainaUe  using  purely  magnetic  focusing.  A  very  simple  lens, 
and  probably  the  optimum  for  this  example,  is  that  of  a  uniform  axi^  field  B  (that  extends 
to  the  image  plane)  for  which  the  chromatic  aberration  coefficient,  Cc,  is  minimized  by 
choosing  the  largest  practical  value  of  B,  about  IT  for  a  conventional  magnet. 


It  is  easily  shown  that,  for  large  demagnification,  the  focal  length  is  given  by 


f  =  -l/k  = 

and  the  chromatic  aberration  coefficient  is  given  by 


For  a  low  voltage  column  the  dominant  aberration  is  chromatic  giving  rise  to  a  disk  of 
confusion  of  diameter: 

dc  =  =-^  “ 

The  diameter  of  the  Airy  disk  (diffraction)  is  given  by: 


The  combined  diameter  can  be  expressed  as: 


d=Vdf2+dc2 


for  3V=leV  and  B=  IT,  aopt=15mrad. 
so  that  dmin  =  I'^dcmin  =  1.4  dfmjii  =  ^  y  ^ 

Thus,  for  purely  magnetic  focusing  with  conventional  magnets  our  lower  limit  of  beam 
diameter  for  V=2000V  is  about  3  nm.  This  is  quite  small  enough  for  our  envisaged  mask 
making.  How  close  we  can  approach  this  limit  when  space  charge  effects  and  finite 
brightness  have  to  be  taken  into  account  is  the  subject  of  current  work. 

Another  configuration  that  has  been  described  for  low  voltage  electron  optics  is  the 
retarding  field  operation.  In  its  simplest  form  a  uniform  retarding  field  is  set  up  between  the 
(flat)  target  and  the  flat  surface  of  the  final  electrode  of  the  beam  forming  system.  If  the 
target  potential,  V  «  electrode  potential,  Ve,  (relative  to  the  source)  then  die  field  is  Ve/d 
where  d  is  the  electrode-to-target  distance  and  most  of  the  focusing  takes  place  in  tfeh 
retarding  field.  Assuming  parabolic  trajectories  the  chromatic  aberration  of  Ais  retarding 
field  lens  is  given  by 

Cc=  d^ 

The  minmum  value  of  d  is  set  by  the  maximum  practical  field,  about  lO^V/m.  Thus  at  best: 

Cc=  V10*7m 

and  dc  =  9V10*^a  m. 


Proceeding  as  before:  ttopt  = 

and  so  the  minimum  value  of  d  is  given  by; 


Thus  for  V=20()0V  and  8V=leV,  dmin  =  2.2  nm  which  is  marginally  better  than  the  result 
with  magnetic  only. 

The  best  result  might  be  expected  to  be  when  both  focusing  fields  overlap.  Again  that  case 
is  being  studied  under  the  current  contract. 

However  it  is  abundantly  clear  that  the  possibility  now  exists  for  obtaining  very  fine  probes 
for  low  voltage  lithography. 
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8.SILIC0N  CALIBRATION  STANDARD  FOR  HIGH  PRECISION 
PATTERN  GENERATORS 


Achieving  accurate  overlay  in  microlithographic  processing  is  one  of  the  key 
requirements  in  pushing  to  higher  density  integrated  circuits.  For  X-ray  and  lx 
lithography,  minimizing  the  errors  is  hampered  by  the  lack  of  a  suitable  reference  against 
which  to  check  mask  and  electron  beam  tool  distortion.  Initial  studies  show  the 
feasibility  of  making  a  high  accuracy  distortion  standard  based  on  the  crystal  lattice  of 
unstressed  Si  wafers.  The  aim  of  the  standard  is  to  facilitate  the  measurement  of 
distortion  by  providing  a  working  reference  of  high  accuracy. 

A  Si  crystal  has  the  property  that  when  cut  and  polished  along  the  (100)  face,  the 
crystal  can  be  anisoiropically  etched  along  orthogonal  [110]  directions  leaving  four  (111) 
faces.  Under  the  right  conditions  lithographically  exposed  lines  will  self  align  to  the 
crystal  axis  and  on  unstrained  Si  can  give  the  straightness  and  orthogonality  required  for 
a  distortion  standard.  Initial  experiments  were  performed  using  a  Canon  stepper  for 
optical  exposures.  Standard  3  inch  wafers  were  exposed  with  a  series  of  lines  and 
squares.  The  pattern  was  steeped  nine  times,  and  at  each  step  the  pattern  was  rotated 
through-0.2®  to  cover  the  wafer  flat  alignment  error.  A  nitride  mask  is  first  etched  with 
the  lithographically  defined  line.  Then  a  KOH  solution  is  used  to  anisotropically  etch  the 
wafer.  The  Si  etches  the  (100)  planes  preferentially  and  is  stopped  by  the  (1 1 1)  planes. 
This  leaves  a  V  shaped  feature  in  the  wafer  which  extends  along  the  (100)  direction,  and 
in  a  perfect  etch  the  feature  has  no  steps  along  the  (1 1 1)  planes. 

Figure  1  and  2  show  SEM  micrographs  from  exposures  with  different  alignments. 
In  figure  1  the  SEM  is  looking  down  into  V  shaped  grooves  comprising  two  (111)  faces 
oriented  in  the  [110]  direction.  The  dark  narrow  lines  are  the  centers  of  the  grooves 
shown  schematically  in  figure  Ic.  In  figure  la  the  wafer  was  well  aligned  and  no  steps 
are  visible.  In  figure  lb  the  series  of  steps  shows  that  the  lithography  is  not  well  aligned 
to  the  crystal  direction,  and  there  are  many  steps  on  the  (1 1 1)  faces.  Figure  2  shows  the 
effects  of  alignment  of  an  etched  square.  Here  the  four  (111)  faces  can  been  seen. 

One  of  the  critical  steps  in  the  production  of  distortion  standard  is  the 
demonstration  that  long  defect  free  lines  can  be  etched  into  Si  and  one  unknown  is  how 
accurately  the  initial  exposed  lines  need  to  be  aligned  to  the  Si  crystal  axis.  For  this 
demonstration  we  have  made  a  Cr  on  glass  photo  mask  in  the  MBES  for  the  contact 
printing  of  lines  on  a  (100)  Si  wafer.  The  lithographic  pattern  consists  of  multiple  lines 
100mm  in  length  arranged  as  a  fan  at  0.4mradian  intervals  over  a  2^  range.  There  are 
three  groups  of  these  lines  and  also  a  set  of  square  features.  The  initial  exposures  are  very 
encouraging.  In  these  first  exposures  it  was  not  realized  that  the  Cannon  printer  skewed 
the  alignment  but  there  was  sufficient  range  in  the  fan  for  the  outer  lines  to  be  etched 
with  a  very  small  number  of  steps  being  detected  in  the  SEM  per  cm.  The  detection  limit 
for  the  size  of  these  steps  being  approximately  lOOnm.  Figure  3  shows  a  micrograph  of  a 
nearly  aligned  line.  The  micrograph  is  of  of  the  V  shaped  line  in  plane  view.  The  dark 
line  running  along  the  center  of  the  micrograph  being  the  bottom  of  the  V,  the  self 
stopping  point.  The  lighter  features  at  the  edges  of  the  V  is  the  nitride  mask,  which  has 
been  undercut.  There  are  very  few  steps  in  the  line  which  shows  the  concept  of  using 


the  fan  of  lines  to  find  the  crystal  axis  does  indeed  work.  The  jagged  edge  of  the  nitride 
mask  in  figure  3  is  from  the  ETEC-MEBES  photomask,  and  is  due  to  the  aliasing  error 
made  in  writing  a  sloped  line.  It  can  be  seen  that  the  anisotropic  etch  ingnores  these 
aliasing  errors.  The  vertex  here  was  imaged  by  inverting  the  SEM  contrast  and  shows 
this  is  a  very  sharp  feature  and  a  good  high  contrast  target  for  electron  beam  metrology. 


Figure  1.  Aligned  and  naisaligned  etched  lines  in  Si. 


a)  Aligned  misaligned 
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c)  Schematic  section  of  lines  anisotropically  etched  into  (100)  Si. 


Figure  2  Square  features  anisotropically  etched  into  Si  (100) 


Figure  3.  Correction  of  MEBES  aliasing  errors. 


The  edge  roughness  of  the  nitride  mask  tech  mask  ( light  feature)  caused  by 
MEBES  aliasing  error  is  not  present  at  the  vertex  of  an  anisotropically 
etched  line  in  the  Si  [100]  direction. 
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INTRODUCTION 

One  of  the  major  difficulties  of  manufacturing  either  optical  or  x-ray  masks  with 
submicTon  features  is  the  proximity  effect.  Electron  scattering  causes  the  proximity  effect 
by  overexposing  dense  features  and  underexposing  isolated  features.  Most  scatter^ 
electrons  fall  within  a  range  of  one  micron  from  the  edge  of  an  electron  beam  during 
exposure  of  the  substrate.  Therefore,  the  proximity  efriect  is  less  important  for  the  large 
feature  masks  that  are  currently  manufactured.  However,  attempts  to  manufacture  both 
optical  and  x-ray  mask,  th  half-  to  quarter-micron  design  rules  has  required  correction 
for  the  proximity  effect. 

At  Etec  we  have  chosen  the  GHOST^  method  to  compensate  for  the  proximity  effect. 
GHOST  works  by  exposing  the  reverse  image  of  a  primary  pattern  at  about  one-half  the 
dose  of  the  primary  pattern,  using  a  large  spot  beam  but  at  the  same  address.  The  result 
smooths  the  energy  density  profile  deposit^  in  the  resist  along  the  edges  of  the  pattern 
and  equalizes  the  dose  between  features.  Figure  1  illustrates  typical  energy  density 
profiles  of  an  uncorrected  primary  exposure,  the  correction  exposure,  and  the  resultant 
corrected  exposure  where  the  two  exposures  are  superimposed. 

The  goals  of  this  project  were  the  following: 

•  Verify  theoretical  GHOST  parameters  for  both  inter-  and  intraproximity  effect 

•  Compare  the  differences  between  patterns  written  with  GHOST  and  patterns  written 
without  GHOST  on  optical  masks. 

•  Quantify  changes  in  the  proximity  effect  to  features  as  small  as  0.3  micron  on  both 
optical  and  x-ray  substrates  with  isolated  opaque  and  clear  lines  from  chip  to  chip  in 
an  array  written  at  optimum  parameters. 

•  Determine  the  minimum  resolution  obtainable. 

•  Determine  linewidth  uniformity  and  line  edge  roughness. 

•  Transfer  patterns  developed  on  a  resist  film  on  x-ray  trilayer  substrates  to  plated  gold 
features.  Characterize,  showing  that  O.S-pm  lithography  can  be  achieved  with 
linewidth  uniformity  of  less  than  or  equal  to  ±5%  of  the  nominal  linewidth. 

RESULTS 

Verification  of  SAMPLE  Parameters  for  GHOST  on 
Optical  and  X-Ray  Substrates 

GHOST  parameters  for  optical  and  x-ray  trilayer  substrates  were  determined  with 
SAMPLE.  The  parameters  were  for  3(X)0A  of  AZ™  5206  resist  on  SOOA  of  midreflective 
chrome  on  a  quartz  substrate  for  optical  masks.  For  x-ray  masks,  the  substrate  was  a 
trilayer  of  30()0A  of  AZ  5206  resist  on  30(A  of  bright  chrome  on  0.7  to  1.0  pm  of  hand- 
bak^  resist  or  polyimide,  on  a  gold/tanu?’  jrn  plating  base.  The  trilayer  is  illustrated  in 
Figiue  2. 

A  6x6  array  dose/correction  ratio  series  was  written  on  both  substrates  with  the  optimum 
theoretical  values  in  the  center  of  the  array.  For  optical  substrates,  the  array  ranged  from  a 
dose  of  7.5  to  9.5  pC  and  a  correction  ratio  of  0.35  to  0.54  QJQp  for  the  optical  masks 
with  a  0. 1  pm  primary  spot  size  and  a  0.58  pm  correction  spot  size.  For  the  x-ray 
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substrates,  the  array  ranged  from  7.5  to  9.0  Qt/Qp  and  0.33  to  0.55  Qc/^  with  a  0.1 
spot  size  for  the  primary  beam  and  a  0.8  lun  spot  size  for  die  correction  Mam.  The  arrays 
were  generated  by  the  GHOST  Process  Evaluation  Tool  (GPET)  program. 

A  process  was  worked  out  for  each  substrate.  Development  was  done  with  undiluted 
PPD  450,  which  allows  iterative  developments.  End  points  were  chosen  by  measuring 
linewidths  of  4.000  ^m  lines  at  several  development  times  on  different  plates  to 
characterize  all  doses.  With  opdcal  masks  the  initial  development  time  was  32  seconds 
followed  by  one  or  more  5'Second  iterations.  The  optical  arrays  were  then  etched,  coated 
with  lOOA  of  gold,  and  examined  in  an  SEM  for  proximity  effect  The  process  for  the 
x-ray  substrates  varied  little  from  that  of  the  substrates;  it  was  34  seconds,  again 
followed  by  various  numbers  of  iteradons.  The  x-ray  substrates  needed  fewer  iteradons  to 
characterize  all  of  the  chips  in  the  array,  and  there  was  a  much  smaUer  proximity  effect 
on  the  chrome/resist  substrate  than  was  observed  on  the  chrome/quartz  substrate.  The 
development  time  required  varied  over  time  for  both  substrates.  A  turnkey  process  is 
being  developed  with  internal  funding  for  GHOSTed  substrates. 

The  optimum  parameters  were  determined  by  examinadon  of  CDs  of  2.000, 1.000, 0.500, 
0.400,  and  0.300  ^m  isolated  opaque  and  clear  lines,  and  a  gradng  structure  to  deteimine 
intraproximity  effect  and  dagger  structures  to  determine  inteiproximity  effects.  Figures  3 
through  5  illustrate  these  features.  The  approach  was  to  find  where  CD-nominal  of  clear 
and  opaque  lines  in  each  of  these  features  begin  to  increase  beyond  normal  noise 
differences.  (In  cases  where  features  were  over-  or  underdeveloped,  the  absolute  value  of 
CD-nominal  for  clear  and  opaque  features  was  compared  and  proximity  effect  was 
determined  when  these  values  were  not  equal  beyond  normal  noise.) 

The  results  showed  that  proximity  effect  was  fully  corrected  to  0.35  iim  at  9.00  |rC  dose 
at  a  0. 1  p.m  primary  spot  size  and  Qc/Qp  of  0.42  at  a  0.58  pm  correcdon  spot  size  on 
opdcal  masks.  There  was  some  distordon  of  opaque  dagger  lines  when  the  clear  feature 
on  either  side  was  0.2  pm.  This  distordon  was  also  observed  in  some,  but  not  all  cases 
with  0.3  pm  clear  features  on  either  side.  The  distortion  was  not  observed  on  other 
features  as  small  as  0.3  pm.  The  difference  is  probably  explained  by  slight  differences  in 
the  MEBES  column  setups  in  the  size  of  the  primary  and  correcdon  spot  sizes,  and  slight 
differences  in  the  dose.  Based  on  these  results,  we  maintain  that  the  proximity  effect  is 
fully  corrected  to  0.35  pm  for  etched  chrome  on  optical  masks.  The  distortion  discussed 
is  observed  in  Figure  6,  a  micrograph  of  an  opaque  dagger  pattern. 

The  results  for  x-ray  substrates  were  determined  by  SEM  evaluation  as  before,  except  that 
we  evaluated  the  gold-coated  resist  rather  than  etched  chrome.  The  results  showed  that  a 
much  smaller  proximity  effect  is  observed  with  this  substrate.  This  is  probably  due  to 
fewer  backscattered  electrons  from  the  mosdy  carbon  substrate.  Ahnost  any  of  the  chips 
in  the  array  could  have  been  chosen  for  the  parameters  by  our  method  of  evaluation.  We 
chose  8.8  pC  at  a  0.1  pm  spot  with  Qc/Qp  set  at  0.50.  At  this  value,  no  distortions  were 
observed  on  any  of  the  features  to  0.3  pm.  Based  on  these  results,  the  lithography  on 
trilayer  x-ray  substrates  is  fully  corrected  to  0.3  pm  and  may  be  good  to  smaller  features. 
Figure  7  is  a  micrograph  of  the  dagger  pattern  in  developed  resist  on  a  trilayer. 

Qualitative  Comparison  of  Corrected  and  Uncorrected  Patterns 

Patterns  of  isolated  lines  and  spaces,  gradng  patterns,  and  dagger  patterns  in  chrome  are 
compared  in  Figures  8  through  11.  Figure  8  compares  corrected  and  uncorrected  clear  and 
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<^)aque  dagger  patterns  for  0.4  4m  lines.  Hie  {voximity  effect  bectxnes  apparent  in  the 
unconected  pattern  beginning  at  about  0.6  pm  and  degrades  as  the  feature  becomes 
smaller.  The  corrected  feature  shows  only  a  small  proximity  effect  for  the  opaque  line  at 
the  0.2  pm  space. 

Figure  9  compares  two  gratings  of  0.3  pm  features.  The  proximity  effect  is  apparent 
without  GHOST.  The  isolated  clear  lines  are  not  developed  out,  die  isolated  t^iaque  lines 
are  narrower  than  the  lines  and  spaces,  the  isolated  t^iaque  pad  is  almost  etch^  away, 
and  the  clear  pads  are  not  etched.  The  cmrected  version  shows  all  features,  each  of  them 
having  s^proximately  the  correa  dimensions. 

Figure  10  shows  a  reverse  tone  of  Figure  9  at  a  higher  magnificadon.  Similar  effects  are 
observed. 

Figure  1 1  shows  isolated  clear  crosses.  Many  of  these,  at  different  sizes,  were  measured 
both  with  an  SEM  and  with  the  EDGE  measurement  system  in  MEBES®  to  quantitatively 
determine  uniformity  and  the  proximity  effect  in  a  3x3-cm  6x6  array  of  chips  written  at 
the  optimum  dose.  The  isolat^  0.3  pm  unconected  line  does  not  develop  out  except  at 
the  cross  where  there  is  a  larger  dose  due  to  backscttered  electrons. 

Quantitative  Measurement  of  the  Proximity  Effect  Measuring  Isolated  Clear 
and  Opaque  Lines  with  SEM  and  EDGE 

A  3x3-cm  array  of  36  chips  containing  clear  and  opaque  crosses  2.000, 1.000, 0.500, 
0.400,  and  0.300  pm  in  diameter  were  written  at  the  optimum  GHOST  parameters  on 
chrome  on  quartz  substrates,  developed  and  etched,  and  measured  with  SEM  and  EDGE. 
The  results  were  determined  by  taking  micrographs  of  the  crosses  as  shown  in  Figure  1 1 . 
Four  measurenxsnts  were  taken  on  each  of  the  36  crosses,  one  from  each  chip. 
Measurements  were  taken  on  both  clear  and  opaque  features.  EDGE  measured  72  lines, 

71  times  each.  Table  1  presents  the  SEM  data  for  both  GHOSTed  and  unGHOSTed 
plates.  Mean  CDs,  mean  deviation  from  nominal,  and  3  sigma  values  are  presented. 

Table  2  presents  the  EDGE  data  for  a  GHOSTed  plate.  Uncotrected  plates  were  not  read 
with  EDGE. 

Figures  12  and  13  plot  the  data  produced  by  EDGE.  Figure  12  is  a  plot  of  the  mean 
linewidth  obtained  from  EDGE  data  versus  feature  size  for  clear  lines  in  x.  All  the  values 
vary  by  less  than  0.012  pm  from  nominal.  The  proximity  effect  that  would  be  seen  with 
simple  change  in  size  is  only  slightly  apparent  at  the  0.3  pm  line,  on  the  order  of  0.01  to 
0.02  pm  correcting  for  the  overdevelopment  of  the  plate.  The  error  bars  are  the  3  sigma 
variation  of  the  lines.  Figure  13  plots  ^e  data  for  opaque  lines  in  x.  All  but  the  0.3  pm 
lines  are  within  0.013  pm  of  nominal.  The  0.3  pm  line  is  within  0.032  pm  of  nominal. 
Correcting  for  overdevelopment,  the  line  is  within  0.02  pm.  The  difference  between  the 
two  lines  shown  in  Figures  12  and  13  can  be  explained  in  terms  of  the  degree  of 
development  ac  etch.  The  lines  were  overdeveloped  or  overetchcd  by  approximately 
0.01  pm.  The  dropoff  of  the  0.3  pm  lines  may  be  the  beginning  of  a  proximity  effect. 
Figure  14  compares  this  data  with  lines  measured  by  SEM  on  an  uncorrected  mask.  The 
solid  lines  represent  the  uncorrected  features  and  the  dashed  lines  the  corrected  features. 
The  amount  of  proximity  effect  observed  between  the  two  is  evident  in  the  slope  of  the 
solid  uncorrect^  lines  as  the  features  get  smaller.  The  results  show  little  proximity  effect 
for  the  corrected  patterns  to  0.3  pm  relative  to  the  uncorrected  values.  For  unconected 
opaque  lines,  the  proximity  effect  is  -0.04,  -0.065,  and  -0.01  pm  for  0.5, 0.4  and  0.3  pm 
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lines,  respectively.  For  clear  lines  ctxrecting  for  overdevelopment,  the  values  are 
-0.06,  -0.1  and  -^.3  for  the  respective  lines.  Corrected  lines  show  a  dropoff  of  about 
0.01  to  0.02  lun  that  could  be  a  result  of  the  proximity  effect  For  this  reason,  we  set  full 
C(»TCCtion  at  0.35  rather  than  0.3  pm. 

Similar  plots  are  shown  in  Figures  15  through  20  for  resist  features  on  the  trilayer.  The 
values  are  very  close  to  ncnninal  for  all  features.  The  0.3  pm  lines  do  not  begin  to  drop  off 
in  these  plots  as  they  do  with  the  chrome  on  glass  substrates.  Therefore,  the  proximity 
effect  is  fully  ccnrected  on  the  trilayer  to  0.3  pm.  The  difference  when  the  two  are 
ctmipared  is  probably  due  to  fewer  backscattered  electrons  in  the  trilayer  due  to  the 
carbon  based  substrate.  This  is  also  bom  out  in  the  large  choice  noted  befrae  when 
determining  parameters  in  the  dose/ccxrection  series  on  the  trilayer.  The  error  bars  are 
maximum  minus  minimum  values. 

Determination  of  Minimum  Resolution 

Features  as  small  as  0.3  pm  in  proximity  patterns  are  relatively  easy  m  produce  with  the 
current  process  for  chrome  on  quartz.  Ctae  can  process  to  4.000  pm  lines  and  produce 
0.3  pm  features  repeatable  where  Unewidth  varies  by  less  than  ±5%  of  the  nominal 
linewidth.  In  about  20%  of  the  cases.  0.2  pm  features  can  also  be  produced.  Figures  21 
and  22  show  0.2  pm  features  in  resist  Figures  22  and  23  show  features  etched  in  chrome. 
When  linewidths  are  mixed  as  in  Figures  22  and  23, 0.2  pm  CDs  are  much  easier  to 
process.  We  believe  that  further  work  may  inqnove  the  resolution  we  can  consistently 
obtain  to  0.2  lines.  Additionally,  proximity  effect  compensation  may  be  improved  by 
using  smaller  spot  sizes  for  the  primary  beam  that  will  be  possible  with  M^ES  IV. 

Lines  of  0.2  pm  were  easy  to  produce  repeatably  in  the  resist  on  a  trilayer  but  were  not 
evaluated.  A  micrograph  is  shown  in  Figure  24. 

Determination  of  Uniformity  and  Line  Edge  Roughness 

The  3  sigma  values  listed  in  Tables  1  and  2  represent  line  edge  uniformities.  It  can  be 
seen  from  the  data  that  the  uniformities  improve  with  the  GHOST  process.  The  3  sigma 
value  is  on  the  order  of  0.04  pm  for  GHOSTed  fearures,  while  it  is  around  0.07  for  the 
uncorrected  features.  This  is  explained  by  the  general  smoothing  of  the  exposure  along 
the  edges  of  the  lines  by  GHOST.  Improvement  of  uniformity  is  also  observed  with  other 
techniques  that  smooth  the  exposure  at  the  edges,  such  as  sizing  the  data  and  MPP™, 
which  introduces  gray  levels  of  exposure  at  the  edges. 

Line  edge  roughness  (LER)  values  determined  by  EDGE  are  listed  in  Table  3.  The  values 
for  all  lines  vary  from  0.037  to  0.055  pm.  The  average  is  0.041  pm.  These  values  are 
typical  of  the  type  of  line  edge  roughness  seen  in  the  cross  depicted  in  Hgure  11.  This 
value  is  within  spec  for  a  MEBES  HI,  but  is  relatively  high  for  features  in  the  0.3  pm 
range.  There  are  several  ways  to  improve  these  values: 

•  The  first  is  to  achieve  a  primary  spot  size  of  as  close  to  0. 1  pm  as  possible.  Normally, 
spots  sizes  are  larger  than  this  by  10  to  20%.  Comparison  of  machine  data  on  spot  size 
has  implied  that  larger  spot  sizes  result  in  greater  line  edge  roughness.  We  have  not 
confirmed  this  by  experiment. 

•  Second,  improve  the  current  process.  A  slower  process  with  no  iterative 
developments  will  also  improve  LER. 
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•  Third,  sizing  the  data  also  improves  the  LER. 

•  Finally,  ocher  resists  such  as  PBS  show  much  better  LER.  Very  smooch  edges  are 
observed  in  the  micrographs  shown  in  Hguxe  25,  which  were  done  with  GHOSTed 
PBS. 

Evaluation  ot  Simulated  X<Ray  Maaka  Written  on  Trilayara 

Five  masks  simulating  x-ray  substrates  were  written  and  sent  to  Du  Prnit,  Danbury,  for 
gold  plating.  They  were  returned  and  evaluated.  Micrographs  are  shown  in  Hgures  26 
through  28. 

The  evaluation  showed  that  a  process  bias  was  introduced  during  the  processing  of  the 
x-ray  substrates  at  Du  Pont  The  bias  varied  from  0.  IS  to  0.2  pm.  In  addition,  the  bias  was 
different  for  clear  and  opaque  features.  For  example.  Table  4  lists  the  mean  linewidths  at 
various  feature  sizes  along  with  the  bias  and  3  sigma  deviations  from  the  mean  for 
32  measurements. 

Table  4  illustrates  the  bias  and  the  shift  in  the  bias  with  the  size  of  the  clear  features. 
Comparison  of  these  values  with  the  values  of  lines  etched  in  chrome  and  resist  lines  on 
the  x-ray  substrates  indicated  that  the  bias  and  shift  is  introduced  in  the  process  during  the 
gold  plating  process  or  RIE.  It  should  be  noted  that  the  substrates  that  were  plated  were 
not  sized  to  account  for  process-induced  bias.  However,  a  plate  was  included  that  had  a 
dose  series  written  on  it  Patterns  that  were  underdeveloped  showed  much  less  bias  and 
shift  from  nominal.  It  is  believed  that  0.5  pm  resolution  can  be  accomplished  through 
gold  plate  with  the  trilayer  process  as  it  is  now.  Uthogr^hic  quality  in  terms  of 
uniformity,  including  3  sigma  values  and  values  of  less  dian  ±5%  of  the  nominal 
linewidth,  may  not  be  accomplished  with  the  process  as  it  is  now. 
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APPENDIX 

Ts^  1 .  SEM  Data  of  GHOSTed  and  unQHOSTed  Patterns 
i  a.  Mean  Unewidths.  GHOSTed 


Nominal  linmrkith 
(mierona) 

Tons 

(as  vahiaa  In  mierona) 

Claar 

Ooaoua 

AHInaa 

1.0 

1.001 

0.992 

0.997 

0.5 

0.505 

0.499 

0.502 

0.4 

0.405 

0.39S 

0.401 

0.3 

0.300 

0.293 

0.296 

b.  Mean  Unewidths.  unGHOSTed 

Nominal  llnawldth 

Tone 

(mierona) 

(an  valuaa  In  mierona) 

Clear 

Opaque 

AHIbiaa 

1.0 

1.073 

0.993 

1.033 

0.5 

0.504 

0.457 

0.481 

0.4 

0.360 

0.335 

0.347 

0.3* 

0 

0.203 

— 

*  0.3  clear  tines  did  not  develop  out 


c.  Mean  Deviation  from  Nominal  Linewidth,  GHOSTed  lines 


Nominal  llnawldth 
(mierona) 

Tone 

(all  valuaa  In  mierona) 

Clear 

Opaque 

Allllnaa 

1.0 

0.001 

-0.008 

-0.004 

0.5 

0.005 

-0.001 

0.002 

0.4 

0.001 

-0.002 

0.001 

0.3 

0 

-0.007 

-0.004 

d.  Mean  Deviation  from  Nominal  Linewidth,  unGHOSTed  lines 

Nominal  llnawldth 

Tone 

(mierona) 

(all  valuaa  In  mierona) 

Clear 

Opaque 

All  llnaa 

1.0 

0.073 

0.5 

0.004 

-0.043 

0.4 

-0.040 

-0.053 

0.3 

-0.300 

-0.097 

— 
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Table  1 .  SEM  Data  QHOSTed  and  unQHOSTad  Patterns  (continued) 
e.  Plus  or  Minus  Percent  of  Deviation  from  Nominal,  QHOSTed 


Nominal  IfeMMrtdth 
(mierona) 

Tona 

Claar 

Obmuia 
■  ■ 

MNnaa 

1.0 

6.8 

0.4 

0.5 

1.0 

0.2 

0.4 

0.4 

0.3 

0.5 

0.3 

0.3 

0 

2.3 

1.3 

f.  Pius  or  Minus  Percent  of  Deviation  from  Nominal.  uhC^OSTed 


Nominal  llnawldth 
(mierona) 

Tona 

a 

Ciaar 

Opaqua 

AHIinaa 

1.0 

7.3 

0.7 

3.3 

0.5 

0.8 

8.6 

3.8 

0.4 

10.0 

16.3 

13.3 

0.3 

100 

313  ' 

— 

g.  Three  Sigma  Deviations,  QHOSTed 


Nominal  ilnawidth 
(mierona) 

Tona 

(all  valuaa  in  mierona) 

Claar 

Opaqua 

All  llnaa 

1.0 

0.034 

0.051 

0.043 

0.5 

0.033 

0.046 

0.4 

0.029 

0.036 

0.3 

0.029 

0.042 

0.037 

h.  Three  Sigma  Deviations,  unGHOSTed 


Nominal  Ilnawidth 
(mierona) 

Tona 

(all  valuaa  In  mierona) 

Claar 

Opaqua 

All  llnaa 

1.0 

0.132 

0.5 

0.082 

0.4 

0.076 

0.057 

0.077 

0.3 

— 

0.060 

— 
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Table  2.  EDGE  Data.  QHOSTed 
a.  Mean  Unewtdths,  GHOSTed 


Nominal  IhMwIdth 
(microns) 

Tons 

(aN  valuas  In  microns) 

Claar 

AH  Unas 

2.0 

2.002 

1004 

2.003 

1.0 

1.001 

0.992 

0.997 

0.5 

0.512 

0.489 

0.501 

0.4 

0.406 

0.395 

0.401 

0.3 

0.295 

0.269 

0.282 

b.  Average  Deviations 


Nominal  llnawldth 
(microns) 

Tons 

(aH  vahisaln  nNerans) 

Claar 

Onaoiia 

AH  Unas 

2.0 

0.002 

0.004 

1.0 

0.001 

-0.008 

0.5 

0.012 

-0.011 

0.001 

0.4 

0.006 

-0.005 

0.001 

0.3 

-0.005 

-0.031 

-0.018 

c.  Percent  of  Deviation  from  Nominai 


Nominal  llnawldth 
(microns) 

Tons 

Claar 

Opaqua 

All  llnaa 

2.0 

0.44 

0.53 

0.48 

0.5 

3.57 

2.66 

3.10 

0.4 

1.52 

1.77 

1.65 

0.3 

1.75 

6.17 

4.25 

d.  Three  Sigma  DeviatiOi  IS 


Nominal  llrMwIdth 
(microns) 

Tons 

(all  valuas  In  microns) 

Clear 

Opaqua 

All  Unas 

2.0 

0.030 

0.034 

0.5 

0.033 

0.056 

0.4 

0.021 

0.027 

0.025 

0.3 

0.012 

0.032 

0.043 
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Tal>te  3.  Line  Edge  Roughness  Determined  by  EDGE 
Mean  3  sigma  values  for  ail  measurements  on  each  line 


Nominal  Hnawidth 
(mlerona) 

Tona 

(all  valuaa  In  mlorona) 

CImt 

Opaqua 

2.0 

0.046 

0.42 

1.0 

0.042 

0.75 

0.041 

0.5 

0.046 

0.038 

0.4 

0.039 

0.042 

0.3 

Oj042 

0.037  ' 
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X-Ray  Masks  —  Materials 
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10.  Appendix  2 

The  Low  Voltage  Alternative  for  Electron  Beam  Lithography 
Y>H.  Lee  et  al.,  Stanford  University 

An  elastic  cross  section  model  for  monte  carlo  simulations  of  low  energy 
electornscatttering  from  high  atomic  number  substrates 
R.  Browning  et  ai.,  Stanford  University. 

Universal  elastic  scattering  cross  sections  for  electrons  in  the  range  /• 
lOOkeV 

R.  Browning,  Stanford  University 

Silicon-on-quartz  reflective  masks  for  0.2S-pm  microlithography 
Y-H.  Lee  etal.,  STanford  University 
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Low  voltage  alternative  for  electron  beam  lithography 

Y.-H.  Lee,  R.  Browning,  N.  Maluf,  G.  Owen,  and  R.  F.  W.  Pease 

Stanford  Solid  State  Laboratory,  Stanford  Uniuenity.  Stanford,  California  943u5 

(Received  24  June  1992;  accepted  3  August  1992) 

The  current  trend  in  electron  beam  lithography  for  patterning  submicron  features  is  towards  the 
use  of  higher  beam  voltages  (20-100  keV).  Among  the  problems  often  perceived  to  be  associated 
with  the  use  of  low  voltages  are  the  poorer  resolution,  the  lower  brightness,  and  the  greater 
sensitivity  to  electric  and  magnetic  interference.  Both  by  simulation  and  by  experiment  at  2  kV 
it  is  shown;  ( 1 )  features  of  less  than  100  nm  are  clearly  resolved  in  resist  of  about  the  same 
thickness;  (2)  such  features  are  clearly  resolved  in  both  sparse  and  dense  pattern;  (3)  such 
features  in  sparse  and  dense  areas  are  clearly  resolved  over  a  twofold  range  of  exposure  doses; 
(4)  such  delineation  is  largely  independent  of  substrate  material;  (S)  there  is  no  evidence  of 
alternating-current  magnetic  interference;  (6)  the  lower  beam  brightness  at  low  voltages  is 
compensated  by  the  increased  sensitivity  of  resists  to  lower  energy  electrons.  The  remaining 
concerns  about  low  voltage  lithography  are  the  reliability  of  resist  with  an  imaging  layer  less 
than  100  nm  thick  and  the  extent  and  effect  of  charging  of  such  a  resist. 


I.  INTRODUCTION 

In  1960,  Mollenstedt  and  Speidel'  used  electron  beams  for 
high-resolution  patterning  on  membranes  by  modifying 
transmission  electron  microscopes.  Practical  electron  beam 
lithography  (EBL)  started  with  a  Westinghouse  team 
(Wells,  Everhart,  and  Matta)^  in  1965,  then  with  an  IBM 
team  (Broers,  Lean,  and  Hatzakis).^  Most  tools  were  mod¬ 
ified  scanning  electron  microscopes  (SEM)  and  primarily, 
for  historical  reasons,  employed  the  same  voltage  ( 10-2S 
kV)  as  did  SEM.^"  Proximity  effects'^  caused  by  lateral 
scattering  of  electrons  of  energies  above  10  keV  in  targets 
have  been  a  subject  of  study  for  about  1 5  yr. 

Two  major  schools  of  thought  have  emerged  to  deal 
with  proximity  effects:  ( 1 )  applying  complicated  proximity 
effect  correction  schemes  to  current  e-beam 
technology, or  (2)  increasing  the  accelerating  voltages 
to  50-100  so  that  the  higher  energy  electrons  for- 

wardscatter  less  in  the  resist  and  the  backscattered  elec¬ 
trons  emerging  over  a  large  area  ( ~  20  ftm  diam )  are  pro¬ 
viding  a  relatively  constant  dose  background.  A  third 
approach  has  been  investigated,  that  is  to  use  low  voltage 
(<2  kV)  such  that  the  scattering  is  laterally  confined  to  a 
small  fraction  of  a  minimum  linewidth.  In  1967,  Pease” 
first  reported  that  by  using  a  retarding  field,  the  chromatic 
and  spherical  aberration  coefficients  can  be  reduced  at  low 
landing  energies  and  that  there  is  no  corresponding  loss  in 
current  density.  In  1986,  McCord  and  Pease^’  used  ul- 
tralow  energy  electrons  (5-100  eV)  in  a  scanning  tunnel¬ 
ing  microscope  (STM)  to  pattern  thin  electron  beam  re¬ 
sists.  It  is  reported  in  this  article  that  using  a  SEM 
operating  at  2  kV  (no  retarding  field)  to  pattern  PMMA 
on  both  high  and  low  atomic  number  substrates;  the  goal  is 
to  quantify  the  expected  improvements  in  critical  dimen¬ 
sion  control  at  low  energy. 


II.  Prediction  of  electron  scattering  range  by 
Monte  Carlo  simulation 

A  Monte  Carlo  simulation  technique  adopted  to  predict 
the  electron  scattering  was  based  on  a  core  program  by 
Joy”  modified  to  employ  an  elastic  cross  section  for  elec¬ 
tron  energy  range  from  1  to  100  keV.^’  The  results  clearly 
indicate  the  smaller  interaction  area  (in  targets)  as  the 
electron  energy  is  lowered  from  10  to  2  keV  (Fig.  1 ).  The 
program  can  also  be  used  to  show  the  distribution  of  en¬ 
ergy  dissipated  per  unit  area  throughout  the  thickness  of 
thin  resist  film  as  a  function  of  radius  from  the  point  of 
impact  (Fig.  2).  As  might  be  expected  this  distribution  is 
independent  of  substrate  materials  (Fig.  2). 

III.  EXPERIMENT 
A.  Exposure  system 

A  Hitachi  S800  SEM  controlled  with  a  PC-based  pat¬ 
tern  generator  (Raith  Elphy  I)  served  as  the  exposure  sys¬ 
tem.  This  SEM  has  a  cold  field  emitter  with  two  lenses,  one 
condenser  lens,  and  one  objective  lens.  The  exposure  pat¬ 
tern  (see  Fig.  3)  used  for  the  study  of  proximity  effects 
consists  of  areas  of  sparse  and  dense  patterns.  Because  the 
SEM  was  operated  in  a  conventional  mode  (no  retarding 


10  kV  JkV  2kV 


Modeluif  system  Mnra  PMMA/Si 


Fig.  I.  Monte  Carlo  simulation  illustrates  that  ( I )  the  lower  the  electron 
energy  the  smaller  the  interaction  area;  (2)  the  smaller  the  lateral  elec¬ 
tron  scattering  range;  (3)  the  thinner  the  surface  imaging  layer 
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Distaiice(nni) 

Fig  2.  Monie  Carlo  simulation  of  electron  energy  dissipated  in  66  nm 
PMMA  on  Si  at  different  accelerating  voltages,  2.  S,  and  10  kV  It  sug¬ 
gests  that  the  lower  the  energy  of  the  electrons,  the  shorter  the  lateral 
electron  scattenng  range.  For  2  kV,  the  curve  applies  for  both  silicon  and 
Au-coated  silicon  substrates 


field),  small  apertures  (  ~  100  fxm  in  diameter)  were  used 
to  reduce  the  aberration  in  the  electron  column. 
Alternating-current  (ac)  interference  was  reduced  by  us¬ 
ing  a  short  working  distance  (5  mm).  The  exposure  doses 
were  estimated  from  the  beam  currents,  which  were  mea¬ 
sured  through  correction  of  the  in  situ  measured  beam 
currents  monitored  via  a  picoammeter.  The  sample  was 
grounded  through  the  picoammeter.  Because  the  in  situ 
measured  beam  currents  were  lower  than  the  total  beam 
currents  due  to  the  backscattered  electrons,  a  Faraday  cup 
was  mounted  on  the  sample  stage  to  measure  the  total 
electron  beam  current  so  that  the  backscattered  electron 
coefficients  can  be  estimated  and  the  in  situ  beam  currents 
corrected.  The  electron  beam  diameter  (43  nm)  was  mea¬ 
sured  by  the  secondary  electron  signal  from  the  Si  knife 
edge  (obtained  by  orientation  dependent  etching)  and  de¬ 
fined  as  the  distance  between  the  10%-90%  points.  The 
address  size  was  controlled  to  be  43  nm. 


two-address  line  in  a  dense  area 

Fig  3  Proximily  effect  pattern  consists  of  two-address  lines  in  dense  and 
sparse  areas  The  dark  areas  are  the  -;.pcsed  areas  The  white  areas  are 
the  unexposed  areas 


Fig.  4.  SEM  of  developed  resist  image  of  86  nm  lines  in  dense  and  sparse 
areas  on  a  silicon  wafer.  The  negligible  linewidth  diff-rence  in  dense  and 
sparse  areas  demonstrated  very  low  proximity  effect  for  2  keV  electrons 
The  white  area  is  the  exposed  area.  The  dark  area  is  the  unexposed  area. 


B.  Resists,  substrates,  and  processing 

Two  types  of  substrates  were  used,  namely,  bare  silicon 
and  silicon  overcoated  with  300  nm  thermally  evaporated 
gold.  PMMA,  66  nm  thick,  was  spin-casted  on  the  sub¬ 
strates  and  prebaked  at  120  ‘C  for  6  h.  The  exposure  dose 
range  was  10-30  fiC/crr?.  After  exposure,  samples  were 
developed  in  1:2  methyl  isobutyl  ketone  (MIBK):isopro- 
pyl  alcohol (IPA)  for  35  s.  Postbaking  was  performed  at 
120  *C  for  20  min.  The  linewidths  were  estimated  by  view¬ 
ing  the  SEM  at  lOOOOOx  (at  25  kV),  although  the 
buildup  of  conumination  in  the  SEM  made  accurate  mea¬ 
surements  difficult. 

IV.  RESULTS  AND  DISCUSSION 
A.  Low  proximity  effect 

Representative  SEM  of  developed  resist  images  of  line 
and  space  patterns  on  a  silicon  substrate  are  shown  in  Fig. 
4.  Note  that  the  linewidth  difference  in  the  dense  and 
sparse  areas  is  less  than  10  nm.  This  experimental  result 
shows  that  proximity  effect  is  very  low  for  2  keV  electrons. 


B.  High  dose  tolerance  and  independence  of 
substrate  materials 

Representative  micrographs  of  the  developed  resist  pat¬ 
terns  on  both  types  of  substrates  are  shown  in  Figs.  5  and 
6.  Note  that  the  difference  of  linewidth  of  two-address  lines 
for  both  dense  and  sparse  patterns  on  both  types  of  sub¬ 
strate  is  negligible.  This  demonstrates  the  virtual  elimina¬ 
tion  of  proximity  effect  on  both  high  and  low  atomic  num¬ 
ber  substrates.  The  allowed  dose  ranges  for  resolution  of  86 
nm  features  in  dense  and  sparse  areas  are  16-28  fiC/ern^ 
for  bare  silicon  and  14-27  /zC/cm^  for  gold-coated  sub¬ 
strate,  respectively.  The  granular  structure  of  gold  made 
precise  linewidth  measurements  difficult. 
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line  in  (  spine  area 

line  in  i  dense  area 


line  in  a  sparse  area 
}  28pC7tm- 
line  in  a  dense  area 


l/2m 


FlCi  S  SEM  of  developed  resisi  image  of  nominal  86  lines  in  dense  and 
sparse  area.t  by  a  2  keV  electron  beam  on  66  nm  PMMA/Si  at  the  expo¬ 
sure  doses  of  lb  and  28  /iC/cm’.  Note  the  nominal  86  nm  lines  can  be 
resolved  in  dense  and  sparse  areas  over  lO'Tr  dose  variation 


Fig  7.  Linewidth  measurement  as  a  function  of  exposure  dose  for  the 
nominal  86  nm  features  on  66  nm  PMMA/Si  expoaed  by  2  keV  electrons 
Note  the  absence  of  any  etTeci  of  the  density  on  the  pattern  and  the  low 
sensitivity  of  linewidth  to  dose  ( A  40  nm  for  a  twofold  increase  m  dose  I 


C.  Dote  effect  on  linewidth 

Measured  linewidths  have  been  plotted  for  both  sparse 
and  dense  patterns  as  a  function  of  exposure  for  both  types 
of  substrate  in  Figs.  7  and  8.  The  measured  linewidth  dif¬ 
ference  ai  a  \0^r  exposure  dose  variation  is  6  nm  on  bare 
silicon  and  9  nm  on  gold-coated  substrates. 

For  linewidth  modeling,  the  following  were  basic  as¬ 
sumptions. 

( I )  The  incident  current  density  distribution  is 


>=yoexp( -r^/of)  .  (1) 

where  r  is  the  radial  coordinate,  measured  from  the  center 
of  the  spot,  and  a,  is  the  radius  at  which  the  current  den¬ 
sity  falls  to  1/e  of  its  peak  value.  The  beam  diameter  d  is 
related  to  a,  by  the  relation  (</=  1.66<t,),  thus  a,  is  25.9 
nm. 

( 2 )  For  a  point  spot,  the  distribution  of  the  energy  de¬ 
posited  in  the  resist  by  the  low-energy  electrons  was  com¬ 
puted  with  the  Monte  Carlo  program.  Therefore,  for  a  spot 
with  radius  the  energy  density  distribution  in  the  resist 


line  in  a  sparse  area 
}  MpC/cm* 

line  in  a  dense  area 


line  in  a  sparse  area 
}  27pC/cm- 
line  in  a  dense  area 


Ipm 

Fl(i  6  SEM  of  developed  resisi  image  of  nominal  86  nm  lines  in  bolh 
dense  and  sparse  areas  by  a  2  keV  electron  beam  on  66  nm  PMMA/.tOO 
nm  Au/Si  al  the  exposure  doses  of  14  and  27  ^C/cm"  Note  that  the  high 
dose  tolerance  applied  to  Au-coated  substrate  as  to  bare  silicon  for  2  keV 
electrons  (see  Fig  5)  The  rough  surface  is  due  to  the  granular  structure 
of  gold 


FUi  8  Linewidth  measurement  as  a  function  of  exposure  dose  of  nominal 
86  nm  features  on  66  nm  PMMA/.t(X)  nm  Au/Si  exposed  b\  2  keV 
electrons 
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Table  I.  Comparison  of  measured  and  modeled  linewidth  difference  at 
10%  dose  variation  for  two-address  lines  in  dense  and  sparse  area  on  both 
St  and  300  nm  thick  gold  coaled  Si  substrates.  Ai.lO'^ :  linemdth  differ¬ 
ence  for  a  10%  exposure  vanauon 


Lines  in  sparse  areas 


Lines  in  dense  areas 


AL10% 

measured 

(nm) 


AL10% 

calculated 

(nm) 


AL10%f 

measured 

(nm) 


AZ.10% 

calculated 

(nm) 


8 

6 

8 

0 

9 

10 

is  a  convolution  of  the  incident  current  density  distribution 
and  the  distribution  of  the  energy  deposited  in  the  resist  for 
a  point  spot. 

Linewidth  variation  was  modeled  for  a  10%  change  in 
exposure  dose  for  two-address  lines  in  dense  and  sparse 
areas  for  both  types  of  substrate.  As  shown  in  Table  1,  the 
measured  linewidth  differences  match  with  the  modeleu 
very  well. 


Fig  9  SEM  of  86  nm  gratings  of  (a)  dryctched  Si  delineated  b>  a  2  ke\ 
electron  beam  along  wiih  dry-etching  processing  (SF^  and  C;CIF.  200 
seem;  150  mTorr,  500  W).  and  (b)  dry-eiched  W  delineated  b\  a  2  keV 
electron  beam  along  wiih  dry-etching  processing  (SF^  and  CBrFi  50 
seem;  50  mTorr.  500  W) 

an  attractive  choice  for  EBL  on  1  x  mask  fabrication  with 
tight  control  of  critical  dimensions. 


D.  Application  of  low  voltage  EBL 

It  was  a  concern  that  low  voltage  electrons  may  not  be 
able  to  expose  resists  which  are  sufficiently  thick  for  useful 
pattern  transfer.  To  investigate  this  potential  problem,  low- 
voltage  ( 2  kV )  EBL  were  employed  for  making  both  op¬ 
tical  and  x-ray  masks.  The  patterned  materials  were  100 
nm  thick  amorphous-Si  (a-Si)  [for  deep  ultraviolet 
(DUV)  optical  masks]  and  200  nm  tungsten  (for  x-ray- 
masks).  The  resist  was  66  nm  thick  PMMA  for  silicon  and 
90  nm  for  the  tungsten.  The  developed  PMMA  resist  lay¬ 
ers  served  as  a  very  effective  plasma  etch  mask  for  a-Si  and 
marginally  effective  for  tungsten.  The  etching  selectivities 
were  measured  by  measuring  the  etching  rates  using  ellip- 
sometry  and  step  profiling.  For  silicon  etching,  we  used  a 
combination  of  SF(,  and  C-CIF,  gases  in  a  Drytek  plasma 
etcher  (both  gases  at  200  seem;  150  mTorr,  500  W).  The 
etching  condition  gives  an  etching  selectivity  of  1;5  for 
PMMA  on  Si.  For  W  etching,  an  etching  selectivity  of  1:2 
of  PMMA  on  W  with  a  combination  of  SF,,  and  CFyBr  in 
the  same  etcher  as  for  silicon  etching  (both  gases  at  50 
seem;  65  mTorr,  500  W)  was  obtained.  Figures  9(a)  and 
9(b)  show  the  dry-etched  silicon  and  tungsten  delineated 
by  the  conditions  described  above.  For  etching  chromium, 
Kuan’*'  has  reported  the  use  of  PMMA  as  thin  as  8  nm. 

V.  SUMMARY 

Low  voltage  (2  kV)  EBL  on  both  high  and  low  atomic 
number  substrates  showed  negligible  proximity  effects  for 
features  down  to  86  nm  as  previously  predicted  by  Monte 
Carlo  simulation.  High  resist  sensitivity  and  a  remarkable 
high  dose  tolerance  were  also  demonstrated.  It  was  also 
demonstrated  that  thin  surface  imaging  layers  are  adequate 
for  pattern  transfer  for  both  optical  and  x-ray  mask  mate¬ 
rials.  Spurious  beam  wobble  due  to  ac  magnetic  fields  was 
negligible.  These  results  suggest  that  low  voltage  (2  kV)  is 
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An  elastic  cross  section  model  for  use  with  Monte  Carlo  simulations  of  low 
energy  electron  scattering  from  high  atomic  number  targets 
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The  Center/or  Integrated  Systems.  Stanford  University.  Stanford.  California  94305 
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The  form  of  the  screened  Rutherford  elastic  scattering  cross  section  has  been  modified  to 
approximate  the  Mott  scattering  cross  section  for  low  energy  electrons  interacting  with  high 
atomic  number  targets,  both  for  the  total  and  the  differential  elastic  cross  section.  A  modified 
form  of  the  total  scattering  cross  section  has  been  found  that  fits  the  Mott  scattering  cross  section 
derived  using  a  partial  wave  expansion.  The  fit  is  from  atomic  number  6-92  over  the  energy  range 
1-100  keV.  The  differential  elastic  cross  section  has  been  modeled  using  a  screened  Rutherford 
distribution  plus  an  isotropic  distribution.  The  ratio  of  forward  to  backscattering  found  from  the 
Mott  cross  section  was  used  as  a  fitting  criteria.  The  screened  Rutherford  distribution  is  fitted  to 
the  forward  scattering  half-angle  of  the  Mott  distribution  and  the  size  of  the  isotopic  distribution 
is  then  adjusted  to  give  the  correct  forward  to  backscattering  ratio.  This  differential  form  has  been 
used  in  a  Monte  Carlo  simulation  of  the  backscattering  from  Au.  Both  the  total  and  the 
differential  cross  sections  are  straightforward  to  implement  in  a  Monte  Carlo  simulation  of 
electron/target  scattering. 


I.  INTRODUCTION 

The  simulation  of  electron  scattering  in  a  target  using  Monte 
Carlo  methods  is  an  important  tool  for  the  understanding  of 
proximity  effects  in  electron  beam  lithography,  and  in  quan¬ 
tifying  electron  microscope  data.  For  these  simulations  it  is 
important  to  have  simple  but  accurate  models  of  the  elec¬ 
tron/atom  scattering  cross  sections.  A  commonly  used  elas¬ 
tic  scattering  cross  section  is  the  screened  Rutherford  cross 
section  which  has  a  convenient  analytical  form  and  is 
straightforward  to  implement  in  a  Monte  Carlo  calculation. 
However,  the  screened  Rutherford  cross  section  can  only  be 
validly  used  for  high  energy  electrons  and  low  atomic  num¬ 
ber  target  materials.  Thus  it  is  inaccurate  for  lithography 
and  scanning  electron  microscopy  which  generally  use  1-20 
keV  primary  electron  energies  and  a  wide  variety  of  target 
materials. 

An  alternative  to  the  screened  Rutherford  cross  section  is 
the  partial  wave  expansion  of  the  Mott  scattering  cross  sec¬ 
tion,  which  is  normally  used  in  the  form  of  a  look-up  table 
and  interpolation  between  data  points.  Here  we  discuss  a 
further  alternative  which  is  to  extract  the  trends  from  the 
Mott  scattering  data  and  apply  corrections  to  the  screened 
Rutherford  formulation.  This  approach  promises  to  give  a 
computationally  efficient  implementation  of  the  Mott  scat¬ 
tering  results. 

The  elastic  cross  section  is  used  in  two  ways  within  a 
Monte  Carlo  simulation,  "^he  total  elastic  cross  section  is 
used  to  define  a  mean  free  path  between  scattering  events, 
and  the  differential  elastic  cross  section  determines  the  angu¬ 
lar  probability  of  scattering.  With  the  screened  Rutherford 
cross  section  the  differential  cross  section  can  both  be  ana¬ 
lytically  integrated  over  a  sphere  to  give  the  total  cross  sec¬ 
tion,  and  it  can  be  integrated  over  a  partial  solid  angle  to  give 
a  simple  form  for  the  scattering  angle  determined  by  a  ran¬ 
dom  number.'  However,  in  this  paper  we  have  treated  the 
total  and  the  differential  cross  sections  separately,  as  it  is  not 
necessary  that  the  same  cross  section  be  used  for  both  cases. 


or  that  a  dimensionless  angular  distribution  be  integrable  to 
an  area  for  the  total  cross  section.  First  we  treat  the  total 
cross  section  and  then  the  differential  cross  section  and  show 
simulation  results  for  the  backscattering  of  electrons  from 
Au  to  illustrate  the  approach. 

II.  THE  TOTAL  ELASTIC  SCATTERING  CROSS 
SECTION 

Without  correction  for  relativistic  effects,  the  total 
screened  Rutherford  cross  section,  cr^  for  the  scattering  of 
an  electron  from  a  single  atom  can  be  written  as 

=  5.21X  10  " • - - - cm-,  (1) 

£-  a(l+a) 

where  Z  is  the  atomic  number  of  the  atom,  E  is  the  incident 
electron  energy  in  keV,  and  a  is  a  screening  parameter  de¬ 
fined  as 

^  U  t>7 

a=  3.4x10- - .  (2) 

E 

It  can  be  seen  that  if  the  screening  parameter  a  is  much 
less  than  1,  which  it  is  over  much  of  our  range  of  interest, 
then  the  screened  Rutherford  cross  section  scales  as  £  -  ' 
and  as  Z  ' Here  we  seek  an  equivalent  expression  for  the 
Mott  cross  sections,  a,, ,  in  terms  of  their  dependence  with  £ 
and  Z.  Reimer  and  Lodding-  have  tabulated  results  of  a  par¬ 
tial  wave  calculation  of  Mott  cross  sections,  and  here  we 
have  plotted  their  results  as  a  function  of  the  reduced  elec¬ 
tron  energy  £=£Z'  Before  reduction,  this  data  was 
first  multiplied  by  the  inverse  relativistic  correction  used  by 
Reimer  and  Lodding. " 

In  Fig.  1,  the  total  Mott  elastic  cross  sections  have  been 
plotted  for  U,  Au,  Mo,  Al,  and  C,  using  the  symbols  U,  A, 
M,  4- ,  and  C.  The  trend  of  the  data  falls  into  two  distinct 
regions:  one  for  reduced  energies  >  0. 1  for  which  C7„  is  pro¬ 
portional  to  £  -  'Z  '  in  agreement  with  the  trends  seen  for 
the  screened  Rutherford  cross  section.  In  the  other,  lower 
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Frc  1  Variation  of  the  total  elastic  cross  section  for  electrons  calculated 
using  the  partial  wave  expansion  ( Ref  2 )  with  the  reduced  electron  energy. 
E  Z  '  "  The  symbols  are  U — uranium,  A — gold.  M — Molybdenum, 
+  — aluminum.  C — carbon 


Fig.  2  Variation  of  total  elastic  scattering  cross  sections  with  electron  ener¬ 
gy  Solid  line.  Eq  (3)  for  L),  Mo.  Al,  C.  and  H  Crosses  are  data  points 
calculated  using  the  partial  wave  expansion  (Ref  2).  The  dotted  lines  are 
the  screened  Rutherford  cross  sections  (Ref  1 ) 


energy,  region  the  results  are  more  complex  with  the  curves 
fanning  out  in  different  trajectories.  For  the  higher  atomic 
number  materials  such  as  Au  and  U  the  trend  suggests  an 
E  dependence  over  an  order  of  magnitude  from  1-10 
keV. 

The  observation  that  there  is  a  systematic  dependence  of 
on  £  and  Z  has  been  made  the  basis  for  finding  a  univer¬ 
sal  expression.  Reimer  and  Lodding’  have  published  results 
or  13  elements,  from  atomic  number  6-92,  and  these  were 
used  to  give  a  fit  across  the  Periodic  Table.  However,  it 
should  be  pointed  out  that  due  to  the  variation  with  atomic 
number  of  the  muffin  tin  potential  used,  Reimer  and  Lod- 
ding’s*  results  are  not  monotonic  with  atomic  number, 
and  the  smooth  fit  discussed  below  ignores  these  variations. 
The  largest  deviations  of  Reimer  and  Lodding’s^  data  from 
the  smooth  fit  are  in  the  region  of  the  transition  elements  and 
will  be  larger  for  the  rare  gas  atoms.' 

The  strategy  for  fitting  the  partial  wave  results  of  Reimer 
and  Lodding'  with  a  form  similar  to  the  screened  Ruther¬ 
ford  cross  section  has  been  described  in  detail  elsewhere.^ 
The  initial  fit  across  the  Periodic  Table  used  the  Mott  cross 
sections  at  100  keV  as  the  normalization  point  and  a  scaling 
law  of  (Z  '  ”  -f  0.032Z')  was  found.  The  region  between 
the  inverse  linear  and  the  inverse  root  regions  was  fitted  us¬ 
ing  a  Gaussian  function  rather  than  using  a  polynomial,  and 
the  final  expression  found  was 

(Z'  "  -f-  0.032Z-) 


£r  =  4.7xl0 


(£-b  0,0155Z'  ”£'") 


X- 


(1  -0.02Z"'c^“  ) 


•  cm* 


(3) 


where  «  =  log,,,  8-£  Z  “  ' 

Figure  2  shows  a  comparison  of  Eq.  ( 3 )  with  the  results  of 
Reimer  and  Lodding*  and  the  screened  Rutherford  cross 
section.  The  continuous  line  is  Eq.  (3)  and  for  U  the  fit  is 
better  than  0.5%  over  the  range  1-100  keV.  The  largest  de¬ 
viations  from  Reimer  and  Lodding’s  data  are  in  the  region  of 
the  transition  elements  where  the  errors  are  up  to  30%, 


largely  due  to  departures  from  the  initial  smooth  fit.  The 
screened  Rutherford  cross  section  is  plotted  as  a  broken  line 
in  Fig.  2  and  is  clearly  inaccurate  for  the  larger  atomic 
numbers. 


ill.  DIFFERENTIAL  ELASTIC  CROSS-SECTION 

As  will  be  shown  below,  the  form  of  the  angular  distribu¬ 
tion  of  scattering  makes  a  significant  contribution  to  the  ac¬ 
curacy  of  a  Monte  Carlo  simulation.  It  has  previously  been 
pointed  out  by  Kotera  et  al.'’  that  the  average  scattering  an¬ 
gle  from  the  screened  Rutherford  cross  section  is  much  larg¬ 
er  than  that  for  the  Mott  cross  section.  The  average  scattered 
angle  for  the  Rutherford  cross  section  can  be  calculated  ana¬ 
lytically  from 

Jo^  sin  6  (d(7/d6)d9 

=  -  .  (  '^  ) 

Jn  sin  9-  (da/d9)d9 

giving 

=  ""Vo ■  \l\  +  a  -rra.  ( 5 ) 

Equation  (5)  is  plotted  in  Fig.  3  along  with  data  taken 
from  Reimer  and  Lodding*  for  Au.  It  can  be  seen  that  the 
average  scattering  angle  using  the  screened  Rutherford  cross 
section  is  always  larger  than  the  corresponding  average  an¬ 
gle  for  the  Mott  cross  section  between  the  energies  1-100  eV. 
Equation  (5),  suggests  the  possibility  of  correcting  the 
screened  Rutherford  formula  to  give  a  similar  average  angle 
as  the  Mott  by  decreasing  the  screening  parameter.  The  for¬ 
mula 

a  =  (0.6 -0.0035  El  a^  (6) 

gives  the  result  shown  as  curve  (c)  in  Fig.  3. 

Although  this  approach  seems  reasonable  this  does  not 
produce  a  scattering  cross  section  that  is  similar  to  the  Mott 
cross  section.  Figure  4  shows  the  form  of  the  differential 
cross  sections  for  1  keV  electron  scattering  from  Au.  In  this 
figure  the  cross  sections  have  been  normalized  at  0°.  The 
screened  Rutherford  differential  cross  section.  Fig.  4(a)  is 
much  broader  than  the  Mott,  Fig.  4(b),  and  it  can  be  seen 
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Ftc.  3.  Average  scattering  angle  as  a  function  of  electron  energy,  for  (a) 
screened  Rutherford,  (b)  partial  wave  (Mott)  (Ref  2).  (c)  the  screened 
Rutherford  using  the  screening  parameter  Eq  (6).  (d)  the  screened  Ruth¬ 
erford  plus  isotropic  distribution  ratioed  using  Eq  ( 8 ) 


that  although  the  average  angle  is  reduced  by  using  Eq.  ( 6 ) . 
Fig.  4(c),  because  the  high  angle  scattering  is  reduced,  the 
small  angle  scattering  is  not  like  the  corresponding  Mott 
form,  and  the  ratio  of  the  forward  to  backscattered  probabil¬ 
ities  are  overestimated  by  a  factor  2. 

As  we  wish  to  fit  both  the  forward  and  backscattering 
distributions  to  give  a  realistic  model  an  alternative  is  to  split 
the  Mott  distribution  into  two  parts,  one  part  for  the  small 
angle  forward  scattering,  and  the  other  the  backscattering 
contribution.  The  forms  for  these  two  parts  of  the  distribu¬ 
tion  were  chosen  to  be  the  screened  Rutherford  cross  section 


Fic.  4.  DifTcrential  scattering  cross  section  as  a  function  of  angle  for  1  keV 
electrons  scatienng  from  Au.  (a)  Screened  Rutherford,  (b)  Partial  wave 
(Mott)  (Ref.  2).  (c)  Screened  Rutherford  using  the  screening  parameter 
Eq.  (6)  (d)  I  he  screened  Rutherford  plus  isotropic  distribution  ratioed 
using  Eq.  (8) 
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and  an  isotropic  distribution.  The  reason  for  this  choice  be¬ 
ing  that  both  forms  give  a  simple  formula  for  the  scattering 
angle  dependence  on  a  random  number. 

The  division  of  the  scattering  distribution  into  forward 
and  backscattering  parts  can  be  justified  by  observing  that  in 
most  simulations  there  is  multiple  elastic  scattering  before 
the  electron  comes  to  rest  or  exits  the  target.  For  example  a 
10  keV  electron  backscattered  from  Au,  undergoes  on  aver¬ 
age  ten  elastic  collisions.  Thus  the  differential  scattering  dis¬ 
tribution  is  self-convolved  several  times,  and  while  the  nar¬ 
row  forward  scattering  half-angle  and  the  ratio  of  the 
distributions  are  preserved,  the  broad  features  of  the  back- 
scattering  distribution  are  averaged  out. 

The  screened  Rutherford  distribution  was  first  modified 
by  changing  the  screening  parameter  to  fit  the  forward  scat¬ 
tering  of  the  Mott  distribution  at  the  half-heights.  An  iso¬ 
tropic  distribution  was  then  added  to  give  the  correct 
weighted  scattering.  Using  the  average  scattering  angle  cri¬ 
teria  produces  the  distribution  shown  as  (d)  in  Fig.  4  This 
distribution  is  now  more  similar  to  the  general  trend  of  the 
Mott  cross  section,  but  underestimates  the  ratio  of  the  for¬ 
ward  to  backscattering  contribution.  A  similar  curve  (not 
shown )  IS  produced  using  the  alternative  criteria  that  the 
ratio  of  the  screened  Rutherford  distnbution  to  the  isotopic 
distribution  is  determined  by  the  ratio  of  the  forward 
(  <  90”)  to  backscattering  (  >  90”)  probability.  The  best  fit  of 
the  half-angle  of  the  screened  Rutherford  distribution  to  the 
Mott  half-angle  is  obtained  using 

a  =  5.5xlO-*Z°'‘V£.  (7) 

The  two  backscattering  criteria  give  similar  trends  for  the 
ratio  of  the  screened  Rutherford  to  isotopic  distributions 
over  the  range  1-100  keV.  For  the  average  angle  criteria 

—  =3.9£:°’^  (8) 

^Iwnropic 

and  the  forward/backscattering  ratio  criteria 


=  4.2£".  (9) 

(^I\olropic 

For  Eq.  (9)  the  explicit  expression  for  the  differential  cross 
section  is  then 


=5.21X10-^'^ 

dil 

.J  1  ,  1) 

id -cos  (9 -a)  4.2£" 


(10) 


where  a  is  defined  in  Eq.  ( 7 ) . 

The  effect  of  using  the  different  criteria  can  be  judged 
from  the  results  of  a  Monte  Carlo  simulation  of  the  back- 
scattering  factor  for  Au  between  0.5  and  10  keV.  The  inelas¬ 
tic  cross  section  used  was  the  Bethe  continuous  slowing 
down  approximation  w  ith  the  Rao-Sahib-Wittry*  low  ener¬ 
gy  extension.  In  Fig.  5  the  upper  curve  (a)  is  the  result  of 
using  the  conventional  uncorrected  screened  Rutherford 
cross  section  [Eqs.  ( 1 )  and  (2)  ].  Using  the  empirical  total 
cross  section  [Eq.  (3)  ],  and  the  differential  screened  Ruth¬ 
erford  cross  section  gives  a  result  that  is  intermediate  be¬ 
tween  the  conventional  screened  Rutherford  model  and  the 
experimental  results.  Using  Eq.  (6).  the  modified  screening 
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Fic.  6.  Variation  of  backscaliering  coefficient  at  normal  incidence  with  elec¬ 
tron  energy  for  Au.  (a)  Screened  Rutherford,  (b)  Screened  Rutherford 
using  the  screening  parameter  Eq.  (6).  (c)  Screened  Rutherford  plus  iso¬ 
tropic  distribution  ratioed  using  Eq.  (8).  (d)  Screened  Rutherford  plus 
isotropic  distribution  ratioed  using  Eq  (9)  Experimental  data  points  A 
(Ref  9)  and  A  (Ref  10). 

parameter,  gives  too  low  a  result,  (b).  The  effects  of  the 
different  criteria  on  the  combined  cross-section  (screened 
Rutherford  plus  isotropic)  can  be  seen  in  the  difference  be¬ 
tween  curves  (c)  and  (d)  of  Fig.  5.  The  backscattering  ap¬ 
pears  higher  between  1-10  keV  for  the  average  scattering 
angle  criteria  than  for  the  forward/backscattering  criteria. 
The  experimental  results  in  this  range  arc  bracketed  by  these 
two  curves  in  this  energy  range,  but  the  forward/backscat¬ 
tering  criteria  gives  results  closer  to  experiment  at  higher 
energies  up  to  100  keV. 

The  differential  cross-section  comprising  the  modified 
screened  Rutherford  distribution  plus  an  isotropic  distribu¬ 
tion  can  be  implemented  very  simply  in  the  standard  Monte 
Carlo  program.  The  choice  of  which  part  of  the  distribution 
to  be  used  at  any  scattering  event  can  be  made  by  an  addi¬ 
tional  random  number  weighting  the  relative  areas  of  the 
distributions  given  by  Eq.  (9).  The  cosine  of  the  scattering 
angle  (a),  is  then  given  by  the  conventional  formula'  for  the 


screened  Rutherford  pan,  and  by  cos  6  =  1  —2R,  where  R 
is  a  random  number  between  0-1,  for  the  isotropic  distribu¬ 
tion. 

IV.CONCLUSiONS 

At  energies  above  0. 1  X  Z  ‘  keV  the  total  elastic  cross- 
sections  derived  from  a  panial  wave  expansion  of  Mott  scat¬ 
tering  of  electrons  from  target  atoms  scales  as  £  “  ‘Z '  At 
lower  energies  the  cross  section  scales  approximately  pro¬ 
portionally  to  £  '  °  Using  these  trends  a  fit  can  be  found 
that  covers  the  range  of  atomic  numbers  from  6-92  and  over 
the  range  1-100  keV. 

A  form  for  the  differential  elastic  cross  section  using  a 
screened  Rutherford  distribution  plus  an  isotropic  distribu¬ 
tion  can  be  found  that  produces  the  correct  magnitude  for 
the  forward  to  backscattering  ratio  over  the  energy  range  1- 
100  keV.  This  form  of  the  differential  cross  sect'on  is 
straightforward  to  implement  in  a  Monte  Carlo  simulation 
of  electron/target  scattering,  and  gives  results  comparable 
to  the  experimental  results  for  the  backscattering  factors 
from  Au. 
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For  high  atomic  numbers  the  total  cross  section  for  electron-atom  elastic  scattering 
calculated  using  the  partial-wave  method  is  found  to  fall  into  two  regimes  over  the  range 
1-100  keV.  The  cross  section  in  the  higher  energy  regime  scales  as  E~  and  in  the  lower 
energy  regime  as  An  empirical  equation  has  been  drawn  up  that  describes  these 

trends,  and  can  give  a  universal  expression  for  the  cross  sections  over  the  Periodic  Table  to 
allow  economical  Monte  Carlo  simulation  of  electron  scattering  over  a  wide  range  of 
energies  and  materials. 


In  the  x-ray  and  Auger  microprobes  it  is  becoming 
important  to  have  a  simple  but  accurate  model  of  the  elec¬ 
tron  scattering  cross  sections  to  perform  quantitative  anal¬ 
ysis  using  Monte  Carlo  simulations.  The  most  commonly 
used  elastic  cross  section  for  this  purpose  is  the  screened 
Rutherford  cross  section.  However,  the  screened  Ruther¬ 
ford  cross  section  is  only  valid  over  a  limited  range  of 
electron  energies  and  materials,  and  becomes  inaccurate  in 
the  range  of  interest  to  scanning  electron  microscope 
(SEM)  analyses  which  generally  use  2-30  keV  primary 
electron  energies.  The  alternative,  which  is  the  use  of  the 
more  accurate  method  of  partial-wave  expansion  using  the 
Mott  scattering  formula,  is  too  computationally  intensive 
to  be  presently  useful  in  on-line  simulations.  Here  we  re¬ 
port  a  universal  total  elastic  scattering  cross  section  de¬ 
rived  from  partial-wave  Mott  scattering  data,  which  like 
the  Rutherford  cross  section  is  simple  to  implement,  but  is 
accurate  over  a  wide  range  of  materials  (carbon  to  ura¬ 
nium)  and  electron  energies  (1-100  keV). 

Without  a  correction  for  relativistic  effects,  the 
screened  Rutherford  cross  section  from  a  single  atom 
can  be  written  as' 


0^  =  5.21X10-^' 


F  ad  -t-a) 


cm 


(1) 


where  Z  is  the  atomic  number  of  the  atom,  E  is  the  inci¬ 
dent  electron  energy  in  keV,  and  a  is  a  screening  parameter 
defined  as 


a  =  3.4xlO-^(Z®‘’VE).  (2) 

It  can  be  seen  that  if  the  screening  parameter  a  is  much 
less  than  1,  then  the  screened  Rutherford  cross  section 
scales  as  and  as  Z' Here  we  seek  an  equivalent 
expression  for  the  Mott  cross  sections  in  terms  of  their 
dependence  with  E  and  Z.  Reimer  and  Lodding^  have  tab¬ 
ulated  results  of  a  partial-wave  calculation  of  Mott  cross 
sections,  and  here  we  have  plotted  their  results  as  a  func¬ 
tion  of  reduced  electron  energy  EsE.Z~  '  ”,  Fig.  1.  Be¬ 
fore  reduction  the  data  were  multiplied  by  the  inverse  of 
the  relativistic  correction  used  by  Reimer  and  Lodding.^ 
There  are  two  distinct  regions:  one  for  reduced  ener¬ 
gies  >0.1,  for  which  is  proportional  to  f-'Z'  ”, 
agreeing  with  the  screened  Rutherford  cross  section.  In 


this  region  the  reduced  total  cross  sections  for  U(92), 
Au(79),  Mo(42),  Al(13),  and  C(6)  arc  similar,  and  close 
to  the  figure  for  the  Rutherford  cross  section.  The  other 
region,  for  reduced  energies  <0.1,  is  more  complex  with 
the  curves  for  the  different  elements  fanning  out  in  dif¬ 
ferent  trajectories.  However,  the  results  suggest  &  E~°^ 
dependence  over  an  order  of  magnitude,  from  1  to  10  keV 
in  Au  and  U. 

The  observation  that  there  is  a  systematic  dependence 
of  on  £  and  Z  is  the  basis  for  finding  a  universal  ex¬ 
pression.  The  results  of  Reimer  and  Lodding^  arc  from  1  to 
KX)  keV,  and  for  13  elements  from  C  to  U.  There  arc  some 
difficulties  in  the  data  set  for  fitting  purposes.  First,  due  to 
the  muffin-tin  potential  used^'^  the  results  are  not  mono¬ 
tonic  with  atomic  number.  The  physical  reason  for  this  is 
the  variation  in  atomic  size,  and  is  most  marked  in  the 
rare-gas  atoms.  ^  Second,  there  are  some  unexpected  trends 
and  inconsistencies.  For  example,  the  cross  section  for  C  at 
1  keV  is  larger  than  would  be  expected  from  the  Ruther¬ 
ford  cross  section  and  the  higher  energy  trends 

A  normalization  point  for  each  atomic  number  was 
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FIG.  1.  Log-log  plot  of  the  total  elastic  scattering  cross  section  for  elec¬ 
trons  calculated  using  the  partial-wave  expansion  ( Ref  2 )  vs  the  reduced 
electron  energy  E.Z"'"  The  symbols  are  U — uranium.  A — gold. 
M — molybdenum,  +  — aluminum,  and  C — carbon 
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first  established  by  using  the  total  elastic  cross  section  at 
100  keV  as  a  base  for  the  trends  towards  lower  energy.  The 
Rutherford  cross  section  scales  as  the  power  1.33  with 
atomic  number  but  this  produces  a  value  that  falls  below 
the  Mott  cross  section  at  high  atomic  numbers.  A  smooth 
fit  to  the  Mott  cross  sections  at  100  keV  can  be  found  using 
ihe  expression  Z'  ”  0.032Z^  This  expression  also  com¬ 

pensates  for  the  departure  from  a  simple  atomic  number 
scaling  law  caused  by  the  energy  scaling  terms  discussed 
below. 

For  one  atomic  number  it  was  straightforward  to  find 
the  fitting  constants  in  the  inverse  linear  and  the  inverse 
root  regions  separately.  However,  finding  a  fit  over  the 
entire  energy  range  using  just  E~  '  and  terms  was 

inaccurate,  and  a  polynomial  using  several  partial  root 
terms  was  found  necessary  to  give  the  correct  transition 
between  the  regions  This  polynomial  was  cumbersome, 
and  as  an  alternative  a  Gaussian  weighting  function  was 
used  to  correct  the  region  between  the  inverse  root  and  the 
inverse  linear  regions.  The  function  used  was 

1 

(1  -0.02Z”V-‘'S’ 

where 

u  =  log;o8.£,Z- '  (3) 

This  function  has  two  constants  and  two  dependencies 
with  atomic  number.  These  dependencies  were  arrived  at 
by  observing  in  the  first  place  that  the  deviation  from  in¬ 
verse  linearity  occurs  at  approximately  the  same  total  cross 
section.  As  the  total  cross  section  scales  in  energy  with 
Z*  ”,  then  the  crossover  point  also  scales  this  way.  In  fact, 
there  is  not  quite  enough  range  in  the  data  to  be  definite 
about  this  trend.  Second,  the  knee  in  the  data  is  shallower 
with  decreasing  atomic  number  and  the  Gaussian  correc¬ 
tion  term  need  not  be  as  large  for  small  Z.  A  square  root 
scaling  with  atomic  number  was  used,  and  although  again 
the  data  does  not  extend  quite  far  enough  to  fix  this  de¬ 
pendence,  the  maximum  correction  at  atomic  number  92 
from  the  Gaussian  term  is  only  20%  and  thus  an  exact 
scaling  towards  lower  atomic  numbers  is  not  critical. 

The  final  problem  is  to  find  the  main  dependence  of  the 
inverse  root  multiplier  term  over  the  entire  periodic  table. 
Scaling  in  the  same  way  that  the  knee  is  scaled,  with  a  Z’ 
dependence,  gives  a  maximum  error  of  3%  over  the  range 
1-10  keV  for  Al,  Mo,  and  U.  Fixing  the  constants  using  the 
best  fit  to  the  U  data  gives  the  complete  empirical  equation 

cr=4.7xl0-’* 

(Z' “ -I- 0.032Z^)  1  , 

^  ( £  -I-  0.0 1 55Z"'£^  -' )  ■  ( 1  _  omTP^e  "  ' 

(4) 

A  comparison  of  the  empirical  equation  (4)  with  the 
screened  Rutherford  and  Mott  results  for  U,  Mo,  Al,  and 
C  is  shown  in  Fig.  2.  Also  shown  is  the  comparison  of  Eq. 
(4)  and  the  screened  Rutherford  cross  section  for  H.  The 
empirical  equation  is  shown  as  a  solid  line,  data  from  the 
Reimer  and  Lodding^  tabulation  without  relativistic  cor- 


E 

u 


o 

u 

e 

o 

O 

JS 

UJ 


o 

H 


Elcctron  Energy,  keV 


FIG  2  Log-log  ploi  of  loial  clastic  scattering  cross  sections  vs  electron 
energy  Solid  line.  Eq  (4)  for  U,  Mo.  Al.  C.  and  H  Crosses  are  data 
points  calculated  using  the  partial-wave  expansion  (Ref  2)  The  dotted 
lines  are  the  screened  Rutherford  cross  sections  (Ref.  1 ) 


rection  is  shown  as  crosses,  and  the  dashed  line  is  the 
Rutherford  cross  section  using  the  form  given  by  Eq.  ( 1 ). 

In  comparison  with  the  tabulated  Mott  data,  the  max¬ 
imum  error  using  Eq.  (4)  for  U  from  I  to  100  keV  is  0.5%, 
while  for  Pb  and  Au  the  maximum  error  is  <  5%  over  the 
range  10-100  keV  and  a  maximum  of  3%  below  10  kev.  It 
can  be  seen  from  Fig.  2  that  the  resultant  fit  for  Mo  in  the 
inverse  root  region  is  good.  The  overall  fit  for  Mo  is  similar 
to  Sb  and  Ag  with  a  maximum  error  of  10%  at  100  keV 
and  3%  at  1  keV.  The  result  for  the  midrange  to  lower 
atomic  numbers  is  not  quite  as  good.  The  maximum  errors 
in  the  range  of  atomic  numbers  from  20-40  is  as  high  as 
30%  with  a  typical  error  in  the  10%  range.  The  source  of 
the  larger  errors  is  the  initial  smooth  fit  over  the  nonmono¬ 
tonic  trend  of  the  cross  sections  at  100  keV. 

The  fitting  error  for  Al  and  Si  is  better  than  4%  over 
the  whole  energy  range  while  the  fit  for  C  is  better  than  3% 
for  5  to  100  keV.  The  fit  is  poor  for  C  below  3  keV,  and 
30%  in  error  at  1  keV.  However,  the  trend  of  the  fit  drops 
slowly  below  the  Rutherford  cross  section  towards  lower 
energies,  as  would  be  expected.  For  the  lower  atomic  num¬ 
bers  the  validity  of  some  of  the  scaling  terms  would  only 
been  seen  in  the  low-energy  region.  Using  the  parameters 
that  match  the  Mott  cross  section  data  that  has  only  been 
presented  graphically  by  Ichimura  and  Shimizu*  (replac¬ 
ing  0.0155  with  0.02  in  the  inverse  root  term),  the  fit  for  Al 
is  reasonable  below  1  keV,  and  within  the  disparity  found 
between  using  the  alternative  Thomas-Fermi-Dirac  and 
the  Hartree-Fock  potentials.* 

In  summary,  at  energies  above  0. 1  x  Z'  the  total  elas¬ 
tic  scattering  cross  sections  Cu,  derived  using  a  partial- 
wave  method  scale  as  £~  'Z'  ^  At  lower  energies  is 
approximately  proportional  to  £“°’.  An  empirical  equa¬ 
tion  has  been  drawn  up  that  is  a  reasonable  description  of 
the  total  elastic  scattering  cross  section  over  the  entire 
periodic  table  with  energies  from  1  to  100  keV.  The  equa- 
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tion  is  simple  enough  for  economic  use  within  a  Monte 
Carlo  simulation.  Further  versions  of  this  equation  might 
well  accommodate  a  much  wider  range  of  energies. 
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Silicon  on  quartz  reflective  masks  for  0.25-(im  microlithography 
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One  approach  to  0.2S-/iin  lithography  currently  being  explored  is  the  unity-magnification 
Markle-Dyson  projection  system.  This  system,  operating  at  >1  =  248  nm,  incorporates  a  reflective 
mask  in  its  design.  This  utilizes  the  internal  reflection  at  the  quartz/film  interface.  To  best  meet 
the  requirements  of  high  reflectivity  and  quarter  micron  processing,  amorphous  silicon  on  quartz 
masks  were  chosen.  Completed  masks  were  used  to  print  0. 19-/im  lines  and  spaces  on  the 
prototype  Markle-Dyson  system,  demonstrating  the  feasibility  of  silicon  reflective  masks  for 
lithography. 
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I.  INTRODUCTION 

Integrated  circuit  lithography  is  typically  done  using  Sx 
reduction  optics;  but  as  feature  sizes  decrease,  correspond¬ 
ing  diifraction-limited  reduction  lenses  become  increasingly 
complicated  to  design  and  to  build.  Projection  optics  design 
can  be  greatly  simplified  by  going  to  a  unity-magnification 
system.  One  such  approach  to  0.25-/tm  patterning  uses  a  two 
element  catadioptic  Markle-Dyson  projection  lens,  having 
an  NA  =  0.7  and  operating  at  /I  =  248  nm. '  This  system  is  a 
working  prototype  for  a  scaled  up  production  stepper.  A  key 
feature  of  the  Markle-Dyson  configuration  is  the  1  X  reflec¬ 
tive  mask.  In  this  article  we  describe  the  design,  fabrication, 
and  testing  of  reflective  masks  for  this  system. 

The  reflective  mask  design,  illustrated  in  Fig.  1 ,  consists  of 
a  flat  quartz  (fused  silica)  substrate  with  a  patterned  reflec¬ 
tive  film  on  the  back  surface.  In  operation,  these  n  asks  uti¬ 
lize  the  internal  reflection  at  the  quartz/film  interface.  To 
improve  object  contrast  and  to  diminish  stray  light  transmis¬ 
sion,  the  patterned  reflective  film  is  coated  with  an  absorber 
film:  this  reduces  the  reflectivity  of  the  spaces  between  the 
reflective  pattern  to  about  1  %  .  In  plan  view,  the  object  half¬ 
plane  contains  reflector  and  absorber  structures,  and  the  ad¬ 
jacent  image  half-plane  is  a  clear  quartz  window.  In  the  fol¬ 
lowing  discussions  reflectivity  of  a  film  refers  to  the 
reflectivity  in  quartz  at  the  quartz/film  interface  at  =  248 
nm.  Normal  incidence  is  assumed  unless  otherwise  noted. 

II.  THICK  FILM  REFLECTION 

The  amplitude  reflectivities  of  a  plane  wave  at  the  bound¬ 
ary  between  two  dielectric  media  of  refractive  indices  n^  and 
n,  may  be  calculated  using  the  Fresnel  equations’ : 

n,  cos  d,  —  n,  cos  9, 

r,  = - ^ ^  ,  ( 1 ) 

n,  cos  ^1  -I-  n,  cos 

rt,  cos  6,  —  /i|  cos  6, 

r  =— c - ^ (2) 

n,  cos  y,  -t-  n,  cos 

Here  i  and  1|  denote  electric  polarization  perpendicular  and 
parallel,  respectively,  to  the  plane  of  incidence;  and  0-,  the 
direction  of  wave  propagation  in  media  1  and  2  respectively, 
are  related  by  Snell’s  law: 


sin  9.  n, 

— ^  =  — .  (3) 

sin  Wj  n, 

For  the  reflective  mask  design,  medium  1  is  a  dielectric;  but 
medium  2  is  absorbing,  having  a  complex  refractive  index 
denoted  by  n,  =  w,  -i- yV,.  For  this  case,  Eqs.  ( 1 ),  (2),  and 
(3)  still  hold  if /I;  is  replaced  by  n,.  By  the  Snell  law  relation¬ 
ship,  the  angle  of  reflection  also  becomes  complex,  and  is 
denoted  here  by  O,;  must  be  replaced  by  B,  in  Egs.  ( 1 )  and 
(2).  The  resulting  complex  Fresnel  equations  are  used  to 
calculate  results  described  in  Sec.  IV. 

For  amplitude  reflectivity  written  in  the  form 
r  =  p  exp(  jS),  the  phase  shift  upon  reflection  is  <5.  At  nor¬ 
mal  incidence  <5.  =  —(5,  =  (5.  The  negative  sign  for  the  par¬ 
allel  polarization  arises  from  rotation  of  the  coordinate  sys¬ 
tem  about  the  transverse  axis  upon  reflection. 


The  intensity  reflectivity  R  is  given  by 

■R;  =  l^il^ 

(4) 

R  =  |r,  !-. 

(5) 

This  gives,  at  normal  incidence. 

(6) 

{/7|  -f  a;,  h"’ 


The  above  equation  holds  for  films  sufticienily  thick  that 


Fic.  I  Applicalion  of  reflecine  mask  in  ihc  Markle-D>son  optical  projec¬ 
tion  s>stem  I  Ref  1  ) 
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negligible  intensity  is  transmitted  to  the  back  side.  For  all 
following  discussions,  “reflectivity”  refers  to  intensity  re- 
flectivity. 

Within  a  medium,  attenuation  of  transmitted  intensity  I, , 
relative  to  its  initial  value  can  be  written  as 


where  z  is  the  depth  into  the  medium  and  d 
tion  depth,  which  can  be  written  as 


(7) 

is  the  absorp- 


A 

4itk 


(8) 


Values  of  reported  refractive  index  data  for  aluminum, 
amorphous  silicon  (a-Si),  silicon,  chromium,  and  novolak 
photoresist  are  listed  in  Table  I,  along  with  calculated  reflec¬ 
tivity,  phase  shift  upon  reflection,  and  absorption  depth. 
These  refractive  index  values  are  used  for  all  calculations 
reported  here.  Refractive  index  values  listed  for  novolak 
photoresist  are  estimates  based  on  observed  behavior. 


III.  REFLECTIVITY  VS  THICKNESS 

To  find  the  reflectivity  for  films  thin  enough  to  transmit,  a 
calculation  for  a  stratified  absorbing  medium  is  required.^ 
Results  from  this  calculation  are  shown  in  Fig.  2,  which 
shows  calculated  plots  of  reflectivity  vs  film  thickness  for 
quartz  coated  with  aluminum,  silicon,  amorphous  silicon, 
chromium,  or  photoresist.  For  the  reflective  mask  geometry, 
incident  light  comes  from  within  the  quartz  (n  =  1.5084), 
and  beyond  the  film  is  air  (n  =  1.0).  The  reflectivities  of  the 
four  reflective  materials  approach  the  bulk  value  at  the  film 
thickness  of  30  nm,  corresponding  to  roughly  five  times  d 
(Sec.  II).  Being  able  to  utilize  thinner  films  allows  one  to 
reduce  etch  selectivity  requirements,  thereby  increasing  pro¬ 
cess  latitude  for  film  patterning. 


IV.  OBLIQUE  REFLECTION 

A  calculation  based  on  the  complex  Fresnel  equations  was 
used  to  study  the  reflectivity  for  a  thick  film  at  oblique  inci¬ 
dence.  Results  from  this  calculation  are  presented  in  Fig  3, 


Table  I  Complex  refractive  index  for  various  materials  at  /i  =  248  nm. 
together  with  corresponding  intensity  reflection  coefficients  with  respect  to 
quartz,  phase  shift  upon  reflection,  and  absorption  depth. 


Material 

n 

K 

R 

5(rad) 

Quartz* 

1. 5084 

0.00 

fl-Si' 

I  66 

3.38 

0.533 

-  2.43 

5.84 

Si' 

1.570 

3.565 

0.573 

-  2.45 

5  54 

Al" 

0  IP 

2.94 

0.901 

-2.20 

6.71 

Cr' 

1.36 

1.91 

0.309 

-  2.08 

10.33 

Air 

1.00 

0.00 

0.040 

0.00 

Photoresist 

1.8 

0.1 

0.009 

-  2.84 

197.4 

•Reference  3. 

”  References  4  and  5. 
•  References  4  and  6. 

References  4  and  7. 
'Reference  8. 


Fic.  2.  Calculated  normal  incidence  quartz/film  reflectivity  vs  film  thick¬ 
ness  for  aluminum,  silicon,  amorphous  silicon,  chromium,  and  novolak 
photoresist. 


which  shows  ( a )  reflectivity  and  ( b )  phase  shift  upon  reflec¬ 
tion,’  both  vs  angle  of  incidence  in  the  fused  silica  medium. 
A  typical  partial  coherence  factor  a  of  0.5  in  the  Markle- 
Dyson  projection  system  ( N  A  =  0.7 )  corresponds  to  a  max¬ 
imum  angle  of  0.23  rad(  =  arcsin{<T*NA/1.5084] ).  At  the 
small  angles  used  for  mask  illumination,  the  splitting  of  per¬ 
pendicular  and  parallel  polarizations  is  minimal.  Since  high¬ 
ly  reflective  materials  for  this  mask  application  will  have 
large  values  of  k,  their  phase  shifts  tend  to  have  similar  val¬ 
ues.  Thus  large  phase  differences,  which  would  be  desired  for 
phase  masks,  cannot  be  obtained  between  highly  reflective 
materials  in  the  same  plane. 


V.  OBJECT  CONTRAST 

Object  contrast  A/q  for  the  reflective  mask  depends  on  the 
relative  reflectivities  of  the  reflector  and  absorber  materials: 


Mo  = 


Rr  +  Ra 


where  and  R^  are  reflectivities  of  reflector  and  absorber, 
respectively.  For  typical  optical  imaging  systems,  image 
contrast  is  a  decreasing  function  of  image  spatial  frequency. 
Since  reliable  imaging  requires  having  adequate  contrast  for 
the  recording  medium  (photoresist),  decreasing  the  starting 
object  contrast  will  therefore  decrease  the  resolution  at  the 
image.  A  reasonable  goal  is  to  have  a  starting  object  contrast 
of  at  least  95%. 


VI.  REFLECTOR  MATERIAL  SELECTION 

Chromium,  aluminum,  and  silicon  were  studied  as  possi¬ 
ble  reflector  materials.  Of  these,  chromium  has  low  reflec¬ 
tance  and  would  offer  marginal  contrast  performance  with 
most  resist-type  absorber  materials.  Aluminum  has  high  re¬ 
flectance  but  is  difficult  to  pattern  to  0.25-//m  geometries. 
Silicon  offers  an  appealing  compromise  in  that  it  has  approx¬ 
imately  55%  reflectivity;  it  has  well-behaved  etching  charac¬ 
teristics;  and  it  has  a  5:1  etch  selectivity  compared  with 
PMMA  resist  (using  the  process  described  in  Sec.  VII).  Us¬ 
ing  a  silicon  reflector  and  a  novalak  photoresist  absorber 
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Fio.  3  ( a )  Quariz/film  interface  reRecuvity  vs  incident  angle,  ( b)  phase  shifts  upon  reflection  vs  incident  angle  Perpendicular  and  parallel  polarizations 
plotted  for  Al,  Si,  Cr,  and  photoresist 


with  a  reflectance  of  \  %  gives  an  object  contrast  of  96%. 

Amorphous  silicon  has  better  etching  behavior  for  fine 
patterning  compared  to  that  of  polysilicon.  Figure  4  shows 
that  polysilicon  films  have  a  rougher  etched  texture  than 
identically  etched  amorphous  silicon  (using  the  process  out¬ 
lined  in  Sec.  VII).  Both  films  are  deposited  on  silicon  test 
wafers.  The  polycrystalline  film  shows  grain  delineation 
arising  from  etch  rate  differences  between  grain  and  grain 
boundary  regions  and  between  different  grain  orientations. 


fh> 


rro  4  Nominal  0.25-/im  line  and  space  pattern  in  (a)  amorphous  silicon 
film  and  ;  b)  polysilicon  film  Both  are  films  on  silicon  wafers  and  are  etched 
under  identical  plasma  etching  conditions.  Finer  microstructure  for  amor¬ 
phous  film  gives  smoother  pattern  delineation  Note  nominal  0  13-/tm  line 
and  space  pattern  produced  in  amorphous  silicon  |  top  half  of  4(a )  ] . 


The  finer  microstructure  of  the  amorphous  silicon  film 
etches  more  uniformly  than  does  a  polycrystalline  film. 

The  best  quality  amorphous  silicon  films  used  were  grown 
using  a  PECVD  process. These  films  were  either  intrinsic, 
or  doped  with  phosphorus.  The  effect  of  doping  did  not  noti¬ 
ceably  affect  measured  reflectivity,  which  was  measured  us¬ 
ing  a  Nanometrics  reflectometer,  for  which  measurement 
uncertainty  is  estimated  to  be  ±  5%.  Reflectivity  measure¬ 
ments  of  quartz  wafers  with  various  amorphous  silicon  film 
thicknesses  are  shown  in  Fig.  5.  These  show  a  55%  saturated 
reflectance  at  thicknesses  greater  than  approximately  30  nm. 
Thus,  films  of  30  nm  or  greater  yield  the  expected  reflectance 


Fig  5.  Quartz/film  refleclivity  measurements  of  amorphous  silicon  vs  film 
thickness. 
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calculated  from  tabulated  refractive  index  data.  This  agree¬ 
ment  indicates  that  interfacial  smoothness  and  film  quality 
are  good.  Patterned  amorphous  silicon  on  quartz  reflective 
masks  were  made  with  film  thicknesses  of  30,  SO.  and  100 
nm. 


VII.  REFLECTIVE  MASK  FABRICATION 

A.  Electron-beam  lithography 

A  quartz  mask  substrate  was  coated  with  an  amorphous 
silicon  film  which  can  be  doped  to  obtain  an  electrical  con¬ 
ductivity  sufficient  for  charge  dissipation  during  electron- 
beam  writing.  This  was  then  spin  coated  with  30-  to  50-nm 
thick  PMMA  electron-beam  resist  and  prebaked  at  120‘‘C 
for  6  h.  Next,  pattern  exposure  was  done  at  10  keV  using  a 
modified  MEBES  1  electron-beam  patterning  system.  After 
exposure,  the  PMMA  was  developed  in  a  solution  of  13 
methyl-isobutyl-ketone  (MIBK):  isopropanol  for  50  s  and 
post-baked  at  120  °C  for  30  min.  Since  the  mask  pattern  is 
written  without  proximity  effect  correction,  the  resist  thick¬ 
ness  is  kept  thin  Incorporating  a  proximity  effect  correction 
scheme  would  allow  thicker  resists  to  be  used. 


B.  Plasma  etching 

With  optimized  etching  conditions,  the  developed 
PMMA  resist  layer  serves  as  an  effective  plasma  etch  mask 
for  a  thin  silicon  film.  A  favorable  combination  of  high  etch 
selectivity  and  low  lateral  etching  was  obtained  using  a  com¬ 
bination  of  SF(,  and  C^F^Cl  gases  in  a  Drytek  plasma  etcher 
(both  gases  at  200  seem;  150  mTorr,  500  W).  These  condi¬ 
tions  give  an  etch  selectivity  of  5:1  for  PMMA  on  silicon, 
calibrated  by  etching  PMMA  masked  silicon  wafers.  This  is 
illustrated  in  Fig.  6,  which  shows  0.4-^m-deep  trenches 
etched  in  silicon  using  a  75-nm-thick  PMMA  etch  mask 


C.  Patterning  the  absorber 

The  final  step  m  the  mask  process  is  lo  pattern  the  absorb¬ 
er  layer  over  the  reflective  pattern.  The  absorber  material  is  a 
standard  g-line  photoresist  ( Shipley  1813)  spin-deposited  'o 
l-1.5-/im  thick  and  patterned  using  conventional  optical 
lithography.  The  absorber  layer  also  serves  as  a  protective 
coating  for  the  patterned  reflector  features. 


VIII.  MASK  INSPECTION 

Figure  7  shows  SEM  micrographs  of  a  reflective  mask 
(with  absorber  removed)  taken  at  25  KeV  These  stuivv 
0  25-,  0.38-,  and  0.5/im  pitch  line  and  space  patterns.  The 
patterned  film  is  100-nm-thick  undoped  amorphous  silicon 
on  quartz,  with  20  nm  of  gold  ev  aporated  onto  the  surface  to 
minimize  SEM  charging  effects  After  inspection  the  con¬ 
ductive  gold  layer  can  be  wet  etched  in  a  solution  of  potas¬ 
sium  iodide  and  iodine  without  harming  the  silicon  features. 
High  performance  optical  inspection  methods  such  as  confo- 
cal  microscopv,  solid  immersion  microscopv.  .md  iieai -field 
microscopy  are  also  being  mvestigaied 


FlO  6  400-nm  deep  trenches  plasma  etched  inio  silicon  using  a  75-nm. 
thick  PMMA  etch  mask,  illusiraiing  a  5  1  etch  seleciivity  Sidewall  fillets 
can  be  minimized  bv  operating  below  the  seleciivilv  limii  ol'  PMMA  as 
shown  in  Figs  4  and  7 


.Nominal  U.  1  3)im 
line  and  space 


Nominal  0.13iim 
line/0. 25pm  space 


Nominal  ().25pm 
line  anti  space 
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Flo  8  (a)  0.2!- and  I  hi  II.  I  V-;/ni  line  iiiK)  N[u.t' pjtii'riiv  til  I  U-/iiii-iliicL 
bi-la>cr  phuutrniM.  madr  from  ihe  proj<xicd  iioagc  of  an  aiiuirphoux  sili¬ 
con  refleciive  mask 


IX.  PROJECTION  LITHOGRAPHY 

A  completed  reflective  mask,  shown  without  absorber 
coating  in  Fig.  7,  was  imaged  using  the  Markle-Dyson  sys¬ 
tem  to  produce  0.25-/im  nominal  line  and  space  features  in  a 
0.5-/im-thick  Shipley  SNR248  negative  resist,  and  0. 19-pm 
nominal  line  and  space  features  in  a  1.0-//m-thick  IBM  bi¬ 
layer  positive  photoresist.'  “  Figure  8  shows  projection  im¬ 
ages  of  (a)  0.25- and  (b)  0  19-/jm  line  and  space  patterns  in 
a  1.0-pni-thick  bi-layer  photoresist. 

X.  CONCLUSIONS 

Modeling  and  experimental  measurements  sho.v  that  sili¬ 
con  films  thicker  than  30  nm  are  well  suited  for  the  reflective 
mask  application.  Of  the  various  types  of  silicon  ftlrns  inves¬ 
tigated.  amorphous  silicon  films  offer  the  best  coinbinaiioii 
of  reneciiviiy  and  piocessing  advantages  fur  line  paiicnnng. 


Measured  reflectance  for  films  at  least  30  nm  thick  are  com 
parable  to  calculated  values.  This  shows  that  high  objecf 
contrast  fora  reflective  mask  is  possible  in  practice.  Working 
reflective  masks  with  0.25 /^m  and  finer  feaiurc,s  were  fabri¬ 
cated  and  successfully  imaged  in  projection.  This  clearly 
demonstrates  the  feasibility  of  using  reflective  masks  for 
high  resblutiori  optical  projection  patterning. 
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